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1. On the supposition that astronomical matter may be treated as incompressible and 
homogeneous, a single star rotating freely in stable equilibrium can be spheroidal or 
ellipsoidal, but of no other shape, while Darwin has shown that both components of 
a binary star must be very approximately of the ellipsoidal shape. 

Both for the interpretation of astronomical observations and for the more general 
purposes of cosmogony, it becomes of importance to examine how the sequence of 
figures assumed by an ideal homogeneous mass will be modified by the compressibility 
and non-homogeneity of actual astronomical matter. 

The general mathematical problem of determining the configurations of stable 
equilibrium of the most general compressible mass is one of great complexity, but 
some important simplifications can be introduced by the use of general considerations. 
These are discussed in §§ 2-5 of the present paper. In § 10 we abandon the general 
problem and turn to a detailed study of the configurations possible when the 
compressibility is such that pressure and density are connected by the law 



— ,^y 



p = K p T — cons., 

where k and y are constants, this of course including the important case of a gas in 
convective equilibrium. The results obtained are summarised, and their astronomical 
bearings discussed, in §§ 45-58. The paper is arranged so that these last sections 
contain the main results of the paper in a form which is free from mathematical 
technicalities ; it is hoped that they will prove intelligible to readers who have 
omitted the more mathematical sections. 

2. For a mass of matter of the most general kind, rotating with angular velocity 
w about the axis of z 9 the equations of relative equilibrium are three of the type 

dp 3Q s /-v 

Sx = P dx-> &C " . . • (1) 

where 

Q = Y + io> 2 (x' + f) ■ • (2) 
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Here V is the total gravitational potential, including that of tidal forces from 
neighbouring stars, if any are present. From equation (l) it follows that the surfaces 
p = cons., Q = eons., and p = cons., all coincide, so that the free surface, being a 
surface of constant pressure, must also be a surface of constant density, say p = <r, 
and must also be one of the equipotentials Q = cons., say Q = C. The value of the 
gravitational potential over the surface must accordingly be 

this being the potential of the rotating mass itself and of certain tide-generating 
masses outside. We may suppose the positions and structure of these tide -generating 
masses to be given, so that V 2 V is given at all points outside the surface of the 
rotating mass, while V is given at all points on the surface by the above equation, 
and vanishes at infinity. By a fundamental theorem in potential theory, it follows 
that V is determined uniquely at every point outside the mass in terms of w, and 
the shape of the boundary, whence dY/dn must be determined uniquely in terms of 
these quantities at every point of the boundary. 
If M is the total mass of the rotating body, 



4ttM 



on 



where the integral is taken over the surface of the body, and oil substituting the 
value of dY/dn, this becomes a linear equation, which may be regarded as determining 
C uniquely in terms of M. 

Thus when M, w and the equation of the boundary are given, it appears that V and 
dY/dn are uniquely determined over the boundary. If the value of p at the boundary, 
say <r, is given, and if we also know the law of compressibility at the boundary, 
then p and dp/dn are uniquely determined over the boundary. Thus not one, but 
two, surfaces of constant density are fixed, namely the boundary S and the surface 
just inside it, say S 7 . 

The mass, say M 7 , inside S 7 is in equilibrium under the rotation w and gravitational 
forces which originate from the external tide-generating masses and also from the layer 
of matter between S and S 7 . It now follows from the preceding argument that the 
surface of constant density next inside S 7 , say S 7/ , is also determined. There are 
now three surfaces of constant density fixed, and by a continual repetition of the 
foregoing process, we can fix all such surfaces in turn. # Thus when the external 
boundary is given, and also M, w and the tide-generating masses, which may be 

* This cannot be defended as a piece of rigorous mathematical reasoning, but there can, I think, be no 
doubt that it is true for practical purposes. I have discussed the mathematical complications elsewhere 
(' Roy. Soc. Proc./ A, vol, 98, p. 413). A more formal proof of the theorem will be found in the i Monthly 
Notices of the R.A.S;,' vol. 77, p, 187. 
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regarded as the data of the problem, the interior arrangement of the matter is 
uniquely determined. 

3. The first consequence of the foregoing theorem is that when the boundary of a 
mass is fixed, all the internal vibrations are necessarily stable. For the change from 
stability to instability can only occur through a vibration of zero- frequency, and this 
would require that there should be two contiguous equilibrium arrangements of the 
interior matter, a possibility which is excluded by the result just obtained. 

The only possible configuration for a mass at rest and under no tidal forces will 
clearly be one in which the boundary is spherical and the surfaces of constant density 
are also spherical and concentric with the boundary. # If the mass is set into slow r 
rotation, this system of concentric spheres will give place to a system of concentric 
spheroids. As the rotation increases further, the surfaces of equal density will no 
longer be strictly spheroidal, but it is clear that there must always be a linear series 
of configurations of equilibrium in which the boundary has the shape of a figure of 
revolution. This series of course reduces to the series of Maclaurin spheroids when 
the matter is incompressible. 

Excluding for the present the case in which a- (the density at the boundary), 
vanishes, it can be shown f that there must be an infinite number of points* of 
bifurcation on this series of figures of revolution, these corresponding, to the different 
sectorial harmonics of the figure. The first point of bifurcation corresponds to the 
second sectorial harmonic ; when this is reached the circular cross-section gives place 
to a slightly elliptic cross-section, and this leads to a series of figures having three 
planes of symmetry and three unequal axes, these figures reducing to the Jacobian 
ellipsoids when the matter is homogeneous and incompressible. It is convenient to 
refer to these two series as the series of pseudo-spheroids and pseudo-ellipsoids 
respectively. The first point of bifurcation on the series of pseudo-ellipsoids corre- 
sponds to a third harmonic displacement and leads to a series of pear-shaped figures. 

It seems almost certain, although it has not been rigorously proved, that this last 
series of figures ends by fission into two detached masses revolving round one another 
as in the ordinary binary star formation. Assuming this, we may regard the passing of 
the first point of bifurcation on the series of pseudo-ellipsoids as the beginning of the 
process of fission. Until this stage is reached we have seen that the only possible 
figures of equilibrium for a rotating mass are pseudo-spheroids and pseudo-ellipsoids. 
For an incompressible mass the pear-shaped figures are unstable, so that spheroids 
and ellipsoids are the only possible figures of stable equilibrium. If the pear-shaped 

* The special application of this to the figure of the Earth has been discussed in a separate paper 
(< Roy. Soc. Proc./ A, vol. 98, p. 413). 

t ' Monthly Notices of the K.A.S./ vol. 77, p. 189. In this paper I was mistaken in thinking that 
fr = presents a true exception to the general theory. The general argument there given failed to prove 
the result in the special case of cr = 0, but the result is true nevertheless, as will appear from the present 
paper. 
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figure proves to be unstable also for all compressible masses, then the only possible 
figures of stable equilibrium for a compressible mass will be pseudo-spheroids and 
pseudo-ellipsoids. 

4. When the mass is compressible a complication can occur which, as it happens, 
does not arise in the incompressible problem. The resultant normal force at any 
point on the surface of rotating mass is — dQ/dn, this being the resultant of gravity 
and centrifugal force. In the incompressible problem, dQ/dn does not vanish, except 
in the unstable configurations at the far ends of the spheroidal and ellipsoidal series, but 
this is not necessarily the case in the compressible problem. We shall find that dQ/dn 
can vanish either on the series of pseudo-spheroids or on the series of pseudo-ellipsoids, 
and when this happens matter will necessarily be thrown off at the points at which 
dQ/dn vanishes. The series of figures of equilibrium, may accordingly be abruptly 
terminated at any stage by the vanishing of dQ/dn. 

5. Thus it appears that the series of figures of equilibrium for a compressible mass, 
until the stage at which fission begins, will consist of pseudo-spheroids and pseudo- 
ellipsoids, these series possibly being abruptly terminated by the vanishing of dQ/dn 
at any point. Our problem is to study these series of configurations ; our method 
will be as follows : — 

For an incompressible mass the density at the centre, which we shall denote by p Q , 
is identical with the density at the boundary, which we denote by a-. For a 
compressible mass, Po will be different from <r, and a rough measure of the extent to 
which the density differs from uniformity will be given by a quantity e defined by 



e — ^ — . . 

Po ■ 



(3) 



We know the solution of the problem when e = ; we require to obtain it for all 
values of e. Our method is to adopt the known solution when e = as generating 
solution and to obtain, by what amounts to a method of successive approximations, an 
expansion for Q in powers of e . The success of the method will depend on the extent 
to which the series so obtained is convergent. 

The value of Q obtained in this way will be a function of x 9 y, z, € and of the 
constants which enter into the law of compressibility. When x 9 y, z are small the 
value of Q will be found to be convergent for all values of e, but as x, y, z increase, 
the range of convergence of the series contracts, But in the most important case we 
shall consider, it is found that the series is convergent, even for points furthest 
removed from the origin, for values of e up to some value between two and three. 
The largest value of e which is of physical interest is e = 1 , corresponding to <r = 0, 
and for this value of e? the series converges with considerable rapidity, so that the first 
few terms will give a fair approximation to the truth. 

Our method, then, is to obtain the series of pseudo-spheroids and pseudo-ellipsoids 
as deformations of the already known series of spheroids and ellipsoids, expanding in 
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powers of the parameter e. The boundaries may accordingly be regarded as distorted 
ellipsoids. 

In a previous paper # I showed how to obtain the potential of a homogeneous mass, 
expanded in powers of a parameter which measured its divergence from the ellipsoidal 
shape. As a preliminary to the present investigation, we must examine the 
corresponding problem when the mass is not homogeneous. 

Potential of a Non-homogeneous Distorted Ellipsoid. 

6. We shall assume, as being adequate for the present problem, that it is possible 
to expand the density in the form 



Po—pB—ps — p* 



(4) 



where p is the density at the origin, which is taken to be the point of maximum 
density, and p 2 , p B9 p 4 , ..., are functions of x, y, z of degrees 2, 3, 4 ..., respectively. 
No terms of degree unity occur because, from our choice of origin, the first 
differential coefficients of p vanish at the origin. The value of p 2 is 



^2 ^2 \ 

1 / 2 O P o . ° P \ 

ox ox oy J 



the differential coefficients all being evaluated at the origin. Our choice of origin 
has been such that p 2 is necessarily positive for all values of x, y and z, so that by a 
suitable choice of direction of axes, it must be possible to express p 2 in the form 



p 2 = u--) (!+!+" 

Let us further put 



where e is a numerical quantity which may for the present be left undefined, but will 
ultimately be taken to be the parameter defined by equation (3), and P is a function 
of x, y and z. Then the general value of p, as given by equation (4), becomes 

P = A)""(po — o") ( L i + V2 + l + € -t > o)> • ....... (5) 

\a h c J 

and the boundary, defined by p = a-, has for its equation 



.2 -, .2 «,2 

a 2 h 2 c 



2 ' 71 ^ J2 ^ fci — ■*-• • ' * \ u / 



When e is not too great this is a distorted ellipsoid. In the paper already referred 
to, I showed how to write down the potential of a homogeneous mass whose 

* < Phil. Trans./ A, vol. 215, p. 27. 
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boundary was given by an equation of the type of (6), the solution occurring as a 
series of powers of e. It is now necessary to extend the method so as to be able to 
write down the potential when the internal density is variable and given by 
equation (5). 

As a matter of convenience, of which the advantage will appear later, we shall 
permit P to include terms of degree 2 as well as terms of degree 3 and higher. 
This is equivalent to regarding the mass under discussion as being arrived at by 
distortion from an ellipsoid of density, 

/ 2 2 2\ 

the only restriction on the distortion being that the point of maximum density must 
remain the origin, and that the density at this point must remain equal to p . But 
distortions which are themselves ellipsoidal are not prohibited, as they would be if 
second-degree terms were excluded from P . 

7. Let q be a function of the density p, defined by 

g 2 -^oz^ ? • ( 8 ) 

so that, as we pass from the centre to the boundary, q varies continuously from 
to 1. The surface of constant density p has for its equation 

2 2 2 

?L + {L 2 + Z ? + eJ> = q\ . (9) 

which may be regarded as a distorted ellipsoid of semi-axes qa, qb, qc. The boundary 
now figures as the special case of q = 1. 

Let the potential of a homogeneous mass of unit density bounded by the surface (9) 
be denoted by V (q) when evaluated at a point outside the surface, and by Y { (q) 
when evaluated at a point inside the surface. 

The whole heterogeneous mass of density given by equation (5) may be regarded 
as split up into shells of constant density, and the potential of the whole will be 
found to be given by the formulae : 

V = <rV (1) -f V (q) dp (10) 

J cr 
*P fPo 

V, = trV<(l) + Y i (q)dp+\ Y Q (q)d P> (11.) 

J cr J p 

the first formula giving the potential at a point entirely outside the mass, and the 
second formula giving the potential at a point inside the mass at which the 
density is p. 
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8. Suppose that P is expressed in the form 






x y z 



2 > 7 2' 2 

XI C 



and let P be a function of x, y, z and a new variable m, given by 



P = F 



q 2 x q 2 y q 2 z \ 



q 2 a 2 + jui' q^+ju' gV + /*, 



so that P is the value of P when ^ = 0. Further introduce £', ?/, f ' defined by 



so that 



£' = x/(q 2 a 2 + /m), &c, 

P = F (g 2 r, gY, gT) 



and /and D (an operator) defined by 



/ 



D = 



* 2 f-5^— + -^ l = (gV+ M )f' 2 +...-l, 



gV + M ' q 2 V + n ' q 2 c 2 +/x 

i \ a 1 



i \ a 2 / i 

o + 



i \ a : 



gW tfa' + Jdf + \q 2 b 2 q 2 b 2 + J d n ' 2 VgV qV+J d£ 



n 



Let <f> (q), a function of £', >/, f ', m and g, be given by 



4>(<?) = 



p_ i/DP+ I/I/YdT- -^— 2 (y/YD 3 P+.. 
4^ 2 2 \4' / / 2 2 . 3 2 \4 l/ / 



1 

8 



2 



DP 2 - -J/D 2 P 2 + -L-/ 2 D 3 P 2 - J— /- 3 D 4 P 2 + 



8 



192 



9216 



+ ^-e 3 [D 2 P 3 

192 



, &c. 



(12) 



When /x = 0, D = and P = P , so that <f> (q) reduces to eP , and 



f+tial-lf^ + yL + i. 

J q 2 q 2 [a 2 If c 2 



~t" 6jl q *~" q J • 



Thus the surface of constant density p is the surface /x = in the family of surfaces 



f+<p(q)lq 2 = (13) 
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From the previous investigation^ it now follows that the value of V (q) is 



Vo{g) 



IT 



.00 



,X> 



f+ 



i>(q) \ (fabc dju 

q 2 M<fa?+r)(qW+v)(&+v)}> ' 



■ (14) 



in which the lower limit of integration // is the root of equation (13) at the external 
point x 9 y, z at which the potential is being evaluated. The same formula (14) 
with // put equal to zero will give the value of V t -(g), the potential at an internal 
point. 

To transform these expressions into a form suitable for use in the present problem, 
introduce new variables X, £ y, f to replace //, £ c/ , v/ y f ', these being defined by 



We now have 



X = M /g 2 ? £=q*i' = x/(a 2 + \) 9 &o. (15) 



P = F(£,, ft 



<ff 



9 9 9 

x + JML. + z 



a?+\ b 2 + \ c 2 +\ 



q\ 



D 



3 



2 



(1_ 



i \ a 2 , / i 

T" 



a 2 + \/df \b 2 (f + Xjdr, 



o + 



i \ a : 



,2 



& + x 



/Sf ? 



while formula (14) becomes 



»0O 



Jx' A 



8 » 



8 % « 



16) 



in which A stands for [(a 2 + X) (6 2 + X) (c 2 + X)] i 5 and X' is the root of 



otf 



+ 



y ■ z 



2 



a 2 + X b 2 + X c 2 4~X 



+ 7rrT+ < £(?) =: 2 2 - • • 



(17) 



The same formula (16) in which X' is put equal to zero will of course give the 
value o£Yi(q). 

9. Attacking equations (10) and (ll), we now have 



'Po 



or 



V (q) dp = ( Po -a) V (q) dq 

Jo 



(Po — *) 



Vo(l) 



' 3 ^Yo (g) r 7^ 



so that equation (10) becomes 



Vo = Po^o ( 1 ) — (po — cr) E , 



• • * » a 



(18) 



* i 



Phil. Trans./ A, vol. 215, p. 27. 
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where 



E,- 



V^W. 



(19) 



This expresses that V is equal to the potential of a homogeneous solid of density 
p minus (p^—a) times the potential E , this last term accordingly representing the 
effect of the falling off of density from the maximum value p . 

Similarly equation (ll) becomes 



"V* = p Vi(l) — (p —o)'E i9 



w 



here 



Ev = 



>q 



<r 



dY (q) 



\2dVM 



<i<f dq2+ \/'W ldq2 



(20) 



(21) 



in which the limit of integration is the value of q at the internal point at which the 
potential is being evaluated. 

The value of V (q) 9 as given by formula (16), is a function of q 2 and A', the two 
being connected by equation (17). Thus we have 



<§» 



dV (q) = 3Vo (q) , dV (q) d\' 
dq 2 dq 2 BX' dq 2 ' 



and the last term vanishes since 3V (q)/d\' contains [q 2 f+<t>(q)~] as a factor, 
we have 

JaA dq 2 I A 



Thus 



dq 



M 



(22) 



and the value of dV { (q)/dq 2 is given by the same formula with A' put equal to zero. 



Using these values we find 



n 



E = irabc q 2 



rX 



V 



1 



dq 2 /A. 



dq 2 , 



(23) 



E/ = Tra&c q 2 

Jo Jx' 



' d<p(q) \d\ 
dq 2 J A . 



'ig 2 + wabc 



,oc 



q' 



• 00 



1 



c<j>(q) \d\ 
dq 2 /A. 



d<f. . (24) 



Both values of E may be regarded as given by a double integration in a plane in 
which A and q are rectangular co-ordinates. Let us first consider E . 
In fig, 1, let PQ represent the curve whose equation is 



x' 






V 



-f 



z 



ar + \ b 2 + \ c 2 + A 



+ <p(q) = q 2 (25) 
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Clearly, when q = 0, the value of X is oo 5 while when q — 1, X has some value \ n \ 
which is the root of 



x 



+ 



tf 



+ -? 



■7 



a 2 +\ ¥ + \ c 2 +\ 



+ <p(l) = 1 



• (26) 



The value of E is obtained by integrating over the area shaded in fig. 1 . 
Changing the order of integration, we find 



»GO 



E, 



IT 



abc 



1- 



90 (q )\ 
dq 2 



dec 



dx. 



(27) 



in which the lower limit q is the root of equation (25), while the lower limit \ ft is the 
root of equation (26). 

The value of E^ is obtained by a double integration in the same plane over an area 
such as that shaded in fig. 2, the different directions of shading distinguishing the 





Ym. I, 



V 2 



areas covered by the two separate integrals in equation (24). 
order of integration, we obtain 



E„- = irabc 



^o 



q \ oq / 



Again changing the 



(28) 



• « • • • 



in which the lower limit q is again the root of equation (25). 

This completes the solution of the potential problem; we now attack the main 
problem of determining configurations of equilibrium. 

General Equations of Equilibrium. 

10. We can only find configurations of equilibrium by assuming a definite law to 
connect pressure with density. We shall accordingly assume the relation 



p == K py — cons 



* : 



. . . (29) 
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The equations of equilibrium (l) now assume the form 



icyp 






3as 3a? 



and these have the common integral 



*T 7-1 

— £ —p r 



y-1 



Q + cons. 



• •»••«* 



(30) 



We have already (equations (3) and (5)) assumed the value of p to be 



Pq 










5 » * * » 



8 t 



• • (31) 



and on expanding p y ~~ l by the binomial theorem, equation (30) becomes 



O = cons. +Kyp ( 



7-1 






X 



2 



3 



i(y-2)(y-3) e2^+e 2 P +... 



• « • 



(32) 



Let the figure of density given by equation (31) be supposed to be a figure of 
equilibrium under a rotation oo and tidal forces of potential V T where 

V T = Tl x 2 +T 2 y 2 +T3Z 2 , (33) 

in which of course t^t^-Vt^ ~ 0. We then have at any point inside the mass 



= po [V,. (1) -eEJ + (t^ + t^ + t^) +^w 2 (x 2 + y 2 ), .... (34) 

and on equating the right-hand members of equations (32) and (34) we obtain a 
solution of the problem. 

11. The simplest solution occurs when the variations of density are small, so that 
e is small ; in this case kj is large. We must not simply put e — in equation (32), 
for otherwise would reduce to a pure constant. We accordingly suppose /eye 
to remain finite, and put e = in all remaining terms. Thus the equation to be 
solved is 



2 2 2 \ 

x . y z x 



Poyi(^) + (riX 2 + T 2 y 2 + T S z 2 ) + iw 2 (x 2 -{-y 2 ) = cona-zcy^-^ + l^ + ^j. . (35) 



z 2 
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The boundary reduces to the ellipsoid 



2 2 2 

2 ' 72 ' ^2 *"" X ? 



a 



c 



and the value of V; ( 1 ), which is the potential of this ellipsoid filled with matter of 

unit density, is 

V,- (l) = — 7rabc (J A x 2 + Jb2/ 2 + J c ^— J), (36) 



where the notation is that I have previously used # in which 



J 



\B,.. 



00 



d\ 



o A(a 2 -h\)(b 2 +\)...' * 



. . . . # . [Oil 



On substituting this value for V t -(l) into equation (35) and equating coefficients 
of x 2 , ?/, z 2 > we obtain, as the conditions for equilibrium, 



Ja- 



j b - 



O) 



6 



7rp abc 2irp abc a 



2' 



t * . • • # • • • 



T 2 



00 



6 



■wpvftbc ivp^abc b 2 * 



(38) 



(39) 



eJ n 



C 



T 3- _ 



6 



in which 



7rp abc o 



2? 



• * • 



t *, • • 



e 



Kyepo 



y-2 



vabc 



. . . (40) 



(41) 



By addition of the corresponding sides of the three equations (36) to (38), we 
obtain 



2 



O) 






abc -wp/jibc 



1 1 



( —, + T~o + ~> ) • 



,a 2 If ' c 2 



(42) 



/ 



It is now clear that equations (38) to (40) together with (41) are simply the 
equations which determine a, b, c, the semi-axes of the ellipsoid which is a figure of 
equilibrium for an incompressible mass. We have, however, found that as far as 
this first approximation the value of p is not necessarily constant throughout the 
ellipsoid; it is given by 



P = Po-(po-o-)l-T 2 + Ti + -g)-. • - • .... (43) 

\CC c 



We have further found the relation connecting 6 with the constants k and y. 
Substituting the value for 6 from equation (41) into equation (32), and equating the 



* * Phil. Trans.,' A, vol. 215, p. 50. 
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value of Q so obtained to the value given by equation (34), the general equation of 
equilibrium is obtained in the form 

Po [V t - ( 1 ) - eEJ + (t x x 2 + r 2 y 2 + t 3 2 2 ) + ^o) 2 (x 2 + y 2 ) 



= cons. —7rp abc9 



2 - 2 + eP ) -^ (y - 2) (s ^ + eP 



a 



a 4 



03' 







+ ^(y-2)(y-3) 2^ 2 +eP 



a' 



. (44) 



12. The first approximation has been obtained by putting c = in this equation; 
we now proceed to higher approximations. A second approximation will be obtained 
by omitting, all power of e beyond the first ; a third approximation by not going 
beyond e 2 , and so on. 

On replacing eP by eP + e 2 Q -f e 8 R -f- ..,, we may suppose that the density expanded 
in powers of e is 



p = Po 



x 



,2 



1-e Z^+6-P -fe 2 Q + e s R +... 



<r 



« « . » « \^hiJ) 



so that the boundary, p = <r, is 



a? 



)« 



2^ - -j- g j7 q -j- e v^q + e Jlvq ~r • • • — 1 • 



• • 



* * 



cr 



. (46) 



The general equation of equilibrium will be obtained from equation (44) on 
replacing eP by 

€ JL q ~T" € I^Jq "T € -Lvq "T" • . . * . • . . . . . . . (4/) 



The value of V z -(l) to be used in equation (42) will no longer be that given by 
equation (36) ; let us suppose the whole value to be 

V < (l) + €AV < (l) + ^V t .(l)+..., . (48) 

this being the internal potential of a homogeneous solid of unit density whose 
boundary is determined by equation (46). Similarly, let the whole value of E,be 
supposed expanded in the form 

Finally, let the value of o> 2 in the complete solution be supposed given by 



w 



27rp abc 



n + e An + e 2 Sn + . . (49) 
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Then the complete equation (44) becomes, on dividing throughout by Trp abc, 

~V-[V,(l)+6AV,(l)+6 2 ^V,(l) + ...~-6(E, + eAE,+ ...)] 

iraOG 

+ ^ 2 + T 2 y 2 + Tjl + U +eAn + € 2 S n +„.)(tf+y*) 

7rp aOC 



cons, —d 



X 



/y«2 



2 ± +eP + e 2 Q + e 3 Ko+ ... -ie(y-2) ( 2 ^ + e P + e 2 Q + • • • 

Ct I \ ct 



X' 



+ ^( y -2)(y-3) S^s + eP + . ...... 



a' 



« a 



(50) 



an equation in which each side is equal to Q/7rp abc. 

13. On equating terms which are independent of e, we of course arrive merely at 
equations (38) to (40) ; these determine a, b, c when the data of the problem are 



given. 



On equating terms in e, we obtain 



abc 



7T 



(AV i (l)-E,.) + A»(.-r ! +2/ 2 ) = -6 P -i( y -2)(2^ 

s, Cv 



X 



2\2~ 



. . . (51) 



In evaluating E< : we put e = 0, and so may neglect </> (q). Thus, from 
equation (28), 



Ev = 7ra&c 



JL 
2 



(i-? 4 ) 



A ' 



• (52) 



in which, from equation (25), q is given by 



9 =-! 



* 2 j. y 2 m ^ 



a 2 4- A 6 2 +X c 2 + A' 



• (53) 



Thus, in the notation defined by equation (37), 

E, = -^7rabc{J AA x i + J BB y^J C c^ + ^ BC y 2 z 2 +...-J). 



(54) 



Using this value for E^, it appears that all terms in equation (51) are of degrees 4 
and 2 except P and AV^(l). The value of AV,(l) depends on that of P , and terms 
of degree n in P give rise to terms of degrees n 9 n — 2, ... in AV,(l). It is accord- 
ingly clear that the solution of equation (51) must be such that P consists of terms 
of degrees 4 and 2 ; the value of AV^(l) is then also of degrees 4 and 2, and, on 
equating coefficients, we can determine the coefficients in P . 

Following closely the notation previously used, # let us assume for P the value 



JL a — 



Lx* Mj/_ 4 Ny 2ly 2 z 2 2mz 2 x 2 2nx 2 y 2 2px 2 ^ 2qy 2 2rz* , r x 

a s b s c s 6V cV a'b" a 4 6 4 c 4 ' * l ° ' 



* * Phil. Trans./ A, vol. 215, p. 54. The small quantity ^e 2 of that paper is replaced by e in this 
paper, otherwise the figure of that paper reduces to the present figure on taking a, ft, y and s all zero. 
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so that (cf. § 8) 

P = L£M~M^ + Nf 4 + 2^T + 2mff ^ . . (56) 

As far as terms of first degree in e, we may, from equation (12), take 

^(g)-=e[P-i/DP+ B V 2 im W 

the remaining terms disappearing because P contains no terms of degree higher than 
four in £ *?, £ The whole potential of the solid of unit density, as far as terms in e, 
is from formula (14) 



so that 



AV,(1)--«M [P-i/DP +A /-DV]f. .... (58) 







in which q is put equal to 1. 
Equation (51) now becomes 



/» X' 



& 



.2\2 



= |(2J A ^; 1 +22J BC? /V-J)-J(y-2)0 2?j) +An(a^+y s ). . . (59) 

\ CI / 

Clearly the left-hand member of this equation will be a linear function of L, M, N, . . ., 
while the right-hand member does not involve these coefficients. Thus the values of 
L, M, N, ... will each be the sum of a number of contributions corresponding to the 
different terms on the right of the equation. Let us suppose that 

x q -— • x q ~t~ x o \ *y ' ^ } ? • • ' * * * * • " \^^/ 

and that 

L = L +L (y— 2), &c, . , (61) 

thus separating the contribution from the term in (y— 2) from the remainder. 

14. It will be remembered that eAV t - (l) is the increase of internal potential 
resulting from deforming the surface of an ellipsoid of unit density so that its surface 

is changed from 

^ 3 /a a + y a /6 3 + 2f 3 /c 3 -l =0, . (62) 

into 

x 2 Ja 2 + y 2 /b 2 + z 2 /c 2 -l + eP = 0. (63) 

Consider the special deformation in which P is of the form 

~2 + 7a + "-§), .......... (64) 

a b c J 
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where \fr is any function whatever. Equation (63) can be solved in the form 

x 2 /a 2 + y 2 /b 2 +z 2 /c 2 -I-y=0, (65) 

and on comparison with equation (62) the deformation is seen to consist of the addition 
of a thin homeoidal shell of uniform density. By a well-known theorem, the potential 
of such a shell is constant at all internal points, so that for the particular deformation 
represented by equation (64), we have AV,(l) = a constant. Incidentally, we may 
note that $V± (l) is also constant, and so on to all orders. 

15. Returning to equation (59), it appears that P'^ must satisfy 

,0 ° J7% / ,,- ,2\ 2 



J [P'-i/DP" +T fc/ a D 8 P"]^-0Po"= -Ws 



« 2 



and from the considerations brought forward in the last section, it is clear that the 
solution is 



\tt 1 / V ^ 

so that 



,2\ 2 



1 — 2 \ ^ ~2 ) ? ' • * • (66) 



L" = fa 4 , &c. ; p" = q" = r" - 0. 
We are left with the problem of determining P', which must satisfy 

[F-i/DF+ 1 fe/ a D 8 F]^-eF 

in*' 

= |-(ZJ AA ^ 4 + 22;J B0 ^V»-J) + An(,x 2 + ?/ 2 ) (67) 

Unfortunately there is no simple means of dealing with this equation, and the 
general solution obtained by direct algebraic treatment is so complicated as to convey 
no meaning at all to the mind. Our method will be to consider first an approximate 
solution of a simple form, this having reference only to the fourth degree terms in 
P ; ; we shall then attempt to estimate the amount of error involved, giving detailed 
calculations and precise, solutions for two configurations of special importance. 

16. The approximate solution we shall consider is 

j^ — _.-- J AA , 6 = —•§■-— Jbc? & c - ....... (68) 

This satisfies equation (67) as regards terms of fourth degree except for the integral 
on the left-hand, so that the error of the solution is roughly measured by the value 
of this integral. 

Now from the definition of the integrals J AA , &c, it is clear, that the greater part 
of the value of these integrals arises from contributions from comparatively small 
values of A, so that to a moderate approximation we shall have 

= ... — k, say (69) 



^AA „ 


_ ^BB 


— "cc . 


__ ^BC 


1 


1 


1 


1 


a 4 


& 4 


c 4 


7 2 2 

c 
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If these relations were strictly true, the value of P ' would be 

p «' - "* 5 ( 2 ^ ■ <70) 

so that, from the consideration of § 14, we should have 



[P'-i/DP +g L/ 2 D 2 P]^ = 0. (71) 

Thus the error is of the order of the error of the approximation (69) multiplied by 
the coefficients of the integral (71), which coefficients are found to be small. It may 
be noticed that equations (69) are strictly true in the spherical configuration in 
which a = b = c, and the error increases as a, &, c become more unequal. 

Thus the error in solution (68) is nil in the spherical configuration ; we shall now 
evaluate it exactly in two other configurations : (i.) the point of bifurcation of the 
spheroidal and ellipsoidal series, and (ii.) the point of bifurcation of the ellipsoidal 
and pear-shaped series. 

17. Ellipsoidal Point of Bifurcation.— At the point of bifurcation of the series of 

spheroids with the series of ellipsoids, the semi-axes a, b, c, and the values of n and 6 

are given by 

a = b = 1*1972, c = 0*69766, 

n = ——=0'lS712 9 = 0*47126, 

2irp 

the scale of length, which is entirely at our disposal, being chosen so as to make 
abc = 1. 

Since the configuration under discussion is now spheroidal, we have 

L = M = n, I = m, p = g, 

so that there are only three coefficients, say, L, m and N, to the terms of fourth 
degree, and two coefficients p, r to terms of second degree. 

For any configuration, formula (48) expresses the internal potential of a solid 
of uniform unit density, so that 

V 2 [V i (l)+eAV i (l)+e^V i (l)+...] = -47r. ...".. (72) 

It follows that AV,(l), (5V,(l), &c, are all spherical harmonics. 
In the present problem in which a = 6, these harmonics are also symmetrical about 
the axis of z, and so are functions of x 2 + y 2 and of z only, so that we may assume 

L\ 

= _ AV,(1) f 

wabc 

= ±c n [(x 2 + tf) 2 -Z{x 2 + y 2 )^ . , (73) 
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in which the qiiantities in square brackets are the most general zonal harmonics of 
appropriate type, and the coefficients c lu d u are chosen so as to be identical with 
those used in my previous paper. Either by direct calculation or by comparison with 
the results obtained in this previous paper, we find 



A — T I T 4 T 1 TT \ 

\ a ar j 



+m 



1 
c 



T , ±_ XT i , 

2 X A 3 C "T 2 2 rL A i C i "T 



8c 



4 



NH 



A'C 2 



• • 



. . . (74) 



id, = 2p(j A A--iI 



r 



AA _ j- 
a / c 



6L T 2L T , m/ oT 
a a 4 c 



1 



A 2 C ~" ~2 -*"AC 



JLmAt — T 
4c 2 \c 2 AC 



\ 

AC 2 /) 



. . (75) 



in which the notation is that already defined in formula (37), supplemented by the 
further abbreviations : 



-AB 



*0G 

'" ~ Jo 



\d\ 



o A(a 2 +\)(b 2 +X)...' 



Jul. ati — 



jTi-XJ • « • 



A(« 2 +\)(6 2 + \) ... ' 



. . (76) 



and the further abbreviations G ABit , F ABi< , will be used when required, to denote 
similar integrals having terms X 3 , X 4 in the numerator. 

18. For computation, it is necessary to construct tables of these integrals. The 
table on the next page contains values which are required both here and later in the 
paper. The method of computation has been fully described elsewhere. # 

Using these values, I find in place of equation (74) 



4c 



ii 



0*010947L-0*064325m + Cr()94707N. ..... (77) 



On equating coefficients of x 4 , x 2 z 2 and z* in equation (59) we obtain 



4c, , — 



L 



'"—""s ~ 2JAA—2 (y-~ 2 ) 

(JO 



16Cn— 



ac 



4 



a 



e 



n _ v/ ^ _ 2 J AC g ^y -/_ 2 .2 



(y~ 2 )r2T2 y • • (^ 8 ) 



N_! 



¥c u -0 i 3 = iJoo-i(y-2) 



ac 







-/ 



•* i 
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Table of Integrals. 







( 


[a= 119723, 


c= 0-69766.) 






J. 


I. 


XI. 


Gr. 


F. 


A 


0-51589 











C 


o- 


96821 










AA 


0' 


22938 


0-18712 








AC 


o- 


47781 


0-28334 








cc 


1' 


05115 


0-45659 


-* 






AAA 


0" 


11850 


0-05953 


0-10177 






AAC 


0" 


26244 


0-10164 


0-13766 






ACC 


o- 


60567 


0-18302 


0-19425 






COO 


1 


44620 


0-34725 


0-28757 






A 4 





06589 


0-02406 


0-02506 


0-06585 




A 8 C 





15204 


0-04451 


0-03786 


0-08335 




A 2 C 2 





•36257 


0-08599 


0-05983 


0*10849 




AC 3 





•88791 


0-17349 


0-09857 


0-14629 




C 4 


2 


•22418 


0-36363 


0-17026 


0-20470 




A 5 - 





•03828 


0-01102 


0-00826 


0.01322 


0-0469 


A 4 C 





•09100 


0-02162 


0-01352 


0-01848 


0-0569 


A 8 C 2 





•22240 


0-04381 


0-02321 


0-02656 


0-0704 


A 2 C 8 





•55496 


0-09243 


0-04102 


0-03978 


0-0893 


AC 4 


1 


•41157 


0-20087 


0-07573 


0-06171 


0-1162 


A 6 





•02290 


0-00545 


0-00321 


0-00366 


• 0080 


A*C 





•05570 


0-01120 


0-00555 


0-00557 


0-0105 


A 4 C 2 


0-13876 


0*02345 


0-01025 


0-00852 


0-0143 


A 3 C 3 


0*35131 


0-05136 


0-01881 


0-01407 


0-0198 


A 2 C 4 


0-90490 


0-11455 


0-03668 


0-02315 


0-0285 


A? 




0-00284 


0-00138 


0-00123 


0-0019 


A 6 C 




0-00608 


• 00247 


0-00201 


0-0027 


A 5 C 2 


™. „ 


0-01294 


0-00496 


0-00314 


0-0040 


A 4 C 3 


,. 


0-02948 


0-00904 


0-00586 


0-0057 


A 2 C 4 




0-06675 


0-01887 


0-00964 


0-0093 



Introducing the value (77) for 4c u , these become three linear equations in L, m 
and N, of which I find the solution to be 



L = M = n = r0273(y-2)-r0466 



a 



l = m = 0-3488 ( y -2)-0-2378 ^ (exact solution). (79) 



N = 0-11845 (y-2)-0-06328 



J 



This may be compared with the approximate solution obtained in § 16, which 

2 A 2 
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may be referred to as "Approximation A." Inserting numerical values, this stands 

as follows : — 

L = M = n = 1-0273 (y-2)-l-0273 



l = m= 0'3488 (y--2)-0*2467 
N = 0*11845 (y~2)-0'0626 



> (Approximation A), (80) 



It will be observed that the error in Approximation A is greatest when y = 2, 
when it is of the order of 2 per cent, of the whole value. 

19. Pear- Shaped Point of Bifurcation. — We proceed now to evaluate the exact 
solution for the configuration at which the Jacobian ellipsoid gives place to the 
pear-shaped figure. At this point we have 



a = 1*88583, 



b = 0*814975, c = 0*650659, 



n 



w 



2 



wp 



= 0*141999, 9 = 0*413607 ; 



the lengths again being chosen so that abc — 1. The values of the integrals 
necessary for the evaluation of the potentials have been given in previous papers. # 
The values of I/, M 7 , N', l\ m\ n 1 are now distinct, and are determined by six 
equations of the type (cf equations (78)) 



4<ai-0 (-t) = £ J aa, &c 



,a 



4Ooo-0 



'23 



21 
Vc 4 



> 



— J«n> &C. 



• ' # » » 



(81) 



J 



in which the coefficients c n , c 22 , c 33 , c 12 , c 2B and c 31 are now distinct, the condition that 
the potential shall be harmonic being expressed by three linear relations connecting 
them. The values of these six potential coefficients have been calculated in my 
previous paper; inserting these and solving the six equations (81) I find the set of 
values given in the second column of the following table, the values given by 
"Approximation A" being given in the adjacent column for comparison : — 



i ,.— ,. „■,.„ — — ■„- ■ ~— — - -—-—-■-—■- y— — — -— — - — — — ■ — — - ■""■ — — -— " — —■-«-— 

Coefficient. 


Exact solution. 


Approximation A. 


Difference per cent. 


1/ 

M' 

N' 

V 

m' 
n' 


-15*4353 

- 0-15560 

- 0-04733 

- 0-08468 

- 0*49850 

- 0-95103 


-10-1768 

- 0*15329 

- 0-04677 

- 0-08424 

- 0-62860 

- 1-18103 


-34 

-1 
-1 

+ 26 

+ 24 



* 



*Phil. Trans./ A, vol. 215, p. 61, and vol. 217, p. 21. 
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Clearly the error is quite large, but is concentrated in the coefficients 1/, m\ n f 
which multiply^terms in x. It is remarkable that the cross section of the figure by 
the plane x = is given accurately to within 1 per cent, by this approximation. 

20. It appears that Approximation A will give a solution accurate to within 
2 per cent, for spheroidal figures, but of error varying from 2 to about 20 per 
cent, for ellipsoidal figures. A tolerable approximation to any required ellipsoidal 
figure could perhaps be got by regarding Approximation A as a first approximation, 
and obtaining the error in this by interpolation between the two errors which have 
been accurately estimated. 

21. A second approximate solution, which has an interest other than that of 
accuracy, may be referred to here. We may use the approximate equations (69) 
to simplify the approximate solution (68), and so obtain a still less accurate 
approximation, which we shall call Approximation B. The approximation is 

L ' = "i a4 ' Z'=-^ 2 ,&c.,. .... . (82) 

and the complete value of P becomes 



p - 1 

X n — o 







J „ .2 J2\ 2 



(y-2)-7 fc+h+j)- ■ W 



,a 2 b c 



In this approximation the quantity h is at our choice ; we must select it so as to 
give as good an average value as possible for the approximately equal quantities 
which occur in equations (69). 

Choosing a suitable value of k 9 I find for the coefficients at the ellipsoidal point of 

bifurcation 

L = M = n = 1-0273 (y-2) -07705 ' 

I = m = 0'3488 (y-2) -0*2616 M Approximation B). . . . (84) 
N = 011845 (y-2) -0*0888^ 

Comparing this with the exact solution given in equations (79), it appears that 
the error is as great as 20 per cent, when y = 2. 

The corresponding approximation at the pear-shaped point of bifurcation is almost 
worthless, the error being one of fully 50 per cent. 

The significance of these approximations will appear later. 

22. We turn now to the evaluation of the terms of second degree in x 9 y y z. 
It has already been seen (p. 172) that p" = q" = r" = 0, so that p 9 g, r are identical 
with p\ q\ r f . Or, in other words, p, q, r do not involve (y-2). 
/Let us suppose that in the general ellipsoidal solution (cf. equation (73)) 

[P_i /"DP+ sV/ 2I)2p ]— = fourth degree terms + 4 (d^ + d^ + d^ 2 ). . (85) 



'o A 
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Then on equating coefficients of x 2 , y* and z 2 in equation (59) we obtain 

Ad,-^- = Aw, (86) 

Ad 2 — tr = An, .......... (87) 

4^-^ = (88) 

c 

The value of 4c?! is readily found to be # 

M = 2i9 J AA -^I AA -^I AB -^I AC + 4 ei , (89) 

Ob o c 



where 

~* TT TT -*-^> TT \ 1 / ^ XX < ^^ J XX 

- Xl A s + — XX AB 2 H - Xl A C 2 I "~*2 I 7F~2 -"-ABC "T 2~2 ^AAC "I §1* 

jot 6 c / \6 c car a b 



^ 2* -1a 8 + 72 -1a 2 B + "T -1-A 2 C (90) 

a b c 

The values of 4c? 2 and 4c? 3 can of course be written down from symmetry. 

It is immaterial whether we insert the values of L, M, N, ... , or of I/, M 7 , N', ... , 
in these equations in evaluating 4e l9 &c, since the contributions from L", M", N", ... , 
which are proportional to a 4 , 6 4 , c 4 , ... , must vanish by § 14. 

23. If we are content to use Approximation B, the solution of our problem assumes 

a very simple form. For I/, M 7 , N', ... are now proportional to a 4 , 6 4 , c 4 , ..., and 

therefore 

e x = e 2 = e 3 = (Approximation B) (91) 

Thus equations (86), &c, reduce to 

2jpJaa-^tIaa-^[Iab- ^Iao--^t = An > ..... (92) 

and two similar equations. 

The quantity An, which has been introduced for convenience only, is entirely at 
our disposal. Let us take An = 0, then the solution of the equations is 

p = q = r = (Approximation B) (93) 

Thus, corresponding to Approximation B, the solution corresponding to rotation w 
given by w 2 J2irp = n has for density 

P -- P .-^$44)-^(l- iY -^4 + ^ . (M) 

* Cf. 'Phil. Trans./ A, vol. 215, p. 60. 
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while the boundary is 



$44^-^(1-^-4 < 95 > 



and so is still ellipsoidal. So far as this approximation goes, it appears that the 
shape of the ellipsoid, as given by the ratio of the semi-axes, depends only on 
a> 2 /27T/o ; it is the same as for a homogeneous mass of uniform density p > equal to the 
density at the centre of the compressible mass, rotating with the same angular 
velocity. 

24. When this somewhat unsatisfactory approximation is abandoned there is no 
means of procedure except to calculate e l9 e 2 and e 3 directly from equations (89), &c. 

Suppose that instead of being given by equation (33) the tidal potential had been 
given by 

V T + eAV T = (Ti + eArO^+^a + eATaJ^ + ^ + eATa)^. . . . (96) 

Then on equating terms independent of e in the principal equation we should obtain 
the same equation as before, but in place of equations (86) to (88), there would be 
three equations such as 

a* 7rp abc 

or, inserting the value of &d x from equation (89), 

2 i >J AA-JlAA-^lAB-JlAc + 46 1 - ^ 4 = An+~^. . . . (97) 

These equations become identical with equations (92) if we take 

Atx = 4t7rp abce l9 &c, (98) 

and the equations then have the simple solution 

p = q == r = 0, An = 0. 

Thus the exact solution can be found by superposing the fourth degree terms 
already calculated on to an ellipsoid which is a figure of equilibrium under an 
additional tidal potential 

e AV T = 4:7rpQabce{e 1 x 2j re 2 'jf J re^), ........ (99) 

this necessarily being harmonic, since e 1 + e 2 + e 3 = 0. 

25. Let a\ b\ d be the semi-axes of a figure of equilibrium under this additional 
tidal force, these differing from the old quantities a, &, c by small quantities of the 
order of e. The equations determining a\ b\ d are (cf. equations (38) to (41)) 

a/ 2 6 f 
J ' A ~ W&V +4e6l ~ 2^W = ^ (1 ° 0) 
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and two similar equations, accented letters all referring to the new figure of 
equilibrium, and the value of 9' being now given by {of. equation (41)) 

e' = ^^ ? (101) 

so that, by comparison with equation (41), 

0V6V = 6abc (102) 

There will, of course, be a solution for the compressible mass corresponding to each 
solution of these equations, which are virtually equations giving a solution of an 
incompressible mass under tidal forces. 

When the main tidal force disappears (r x = r 2 = t 3 = 0) the solutions will corre- 
spond, except for small tidal terms, to Maclaurin spheroids and Jacobian ellipsoids. 
The compressible solution which corresponds to the Maclaurin-Jacobian point of 
bifurcation will represent a point of bifurcation for the compressible mass. Let us 
proceed to determine this point of bifurcation exactly. 

At the point of bifurcation the configuration is spheroidal, so that e x = e 2 . By 

subtraction of corresponding sides of equation (100) and the similar equation in 6, 

we obtain 

' 1 1 

d 2 V 
or 

(a /s -6 /a )(J , AB -0'/a /1 6 /a, ) = O (103) 



V V _ at 



The spheroidal series is determined by the vanishing of the first factor, and the 
ellipsoidal series by the vanishing of the second factor. At the point of bifurcation 
both factors vanish, so that a 1 = b\ and 

a J AA = 6 , 
or again, using relation (102), 

^c'J'aa = 0a 2 c (104) 

This equation, together with (100) and its companion, will determine the values 
of a\ d and w at the point of bifurcation. 

26. The values of these quantities differ by terms of the order of e from the 
corresponding quantities a, 6, c at the Maclaurin-Jacobian point of bifurcation, 
so that we may suppose that 



a n a 2 \a 



, 2 = ^+eA(^Uc, / = (o 2 + eA(o 2 (105) 



The equation of the figure of equilibrium is now 



^''■-W©^®^... 



2 * T~2 ""•" e 

a b 
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and this agrees exactly with the figure previously considered if we take 






2r 



\c 



(106) 



Using relations (105) and (106), equation (100) transforms to 



abc J A — 



Ja~ 


Tl 


el 


2 -"-AA J " 

. a I 


r j 

2 -*-AC "" 
G 


Aco 2 

27rp abc 




7rp abc 
+ eabc 


27rp abc a 2 _ 
4e 1 +£>(2J AA - 


_2pd 

4 

a 



= 0. . 



(107) 



The first line vanishes in virtue of equation (38) and the equation reduces to 



4e 1 +p(2J AA -~2l AA 
\ a 



r T 2p6 

~~2 ■*■ AC " i~ 

c a 



An. 



* • * 



(108) 



which agrees with equation (86), as of course it ought to do, and there is a similar 
equation agreeing with equation (88). 

27. Equation (104) may be written in the slightly symbolical form 

a 6 cJ AA + eA(a 6 cJ AA ) = 6a 2 c } (109) 

and the condition which determines the ordinary Maclaurin-Jacobian point of 
bifurcation is 

a Jaa = & 

Thus equation (109), which expresses the condition to be satisfied at the point 
of bifurcation on the compressible series, reduces to 

A (a 6 c J AA ) = 0. . , . . . . . . . . (HO) 

Now from the definition of J AA we have 



»0O 



a 6 cJ 



A.A 



d\ 



o(l+A/a 2 ) 3 (l4-A/c 2 )*' 



whence we readily find that equation (110) is equivalent to 



61 A s( —5) +1 A 2 C ( "2 ) — 0. 



-A 3 \ 2 



# * ♦ # 



. . . (Ill) 



Equation- (ill), together with equation (108) and its companion, determine the 
values of p, r and An at the point of bifurcation on the series of compressible 
configurations. 
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28. The material for numerical computations has already been given in § 18. 
Using the exact solution there given, I find, by direct computation from equation (90), 

4^ = 4e 2 = -0*03404 ; 4e 3 = 0*06808 (112) 

From equation (ill), and the companion to (108), I now find 



P- = -0*016037; 

a 



r 

~2 



c 



0*056337, (113) 



and equation (108) now gives 



An = 



Aw 



27rp ahc 



0*04400. . 



« » 



» t 



. . (114) 



These quantities are all comparatively small ; they would all have vanished under 
Approximation B. 

Second Order Solution. 

29. We now pass to the solution of higher order still. On equating terms in e 2 in 
equation (50), we obtain 

—^ (§Y { (1 ) - AE,) + Sn (x 2 + y 2 ) 

TtCOOG 



= -e 



X 2 \ . i / „s / „\L^ W 



Q -( y -2)P 2^- +i(y-2)(y-3) S±i 



a 



a 



• • 



(115) 



30. From equation (28), taking terms as far as first powers of e only, 



E,- + AE,- = irabc 



J0 [_Jq \ CQ 



dX 

A 



(116) 



in which the lower limit q is given by 



z 2 



+ t™— + -7~~r +j> (q) = q 2 



a 2 ~hX h 2 + X ' o 2 -f-X 



(117) 



As far as e, equation ( 12) gives 



while from § 8. 



^(g) = e[P-£/I)P+^/ 2 D 2 P],. 



• • « • • 



. (118) 



2-/' 



qV 



x 



.2 



+ ^- + ^--q 2 (H9) 



a s +\ b 2 + \ c 2 +\ 



D 

T 



1 1 






i \ a 



• + /i 



i \ a 2 



'"N. £'-> 



a 2 a* + \/dg a \b 2 b 2 +\Jd v 2 \c 2 <f + \/ci 



... . (120) 
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The value of D/g 2 is Independent of g, and $ (q) may be rearranged in the form 



4> (q) = 
whence, on differentiation 



q w ' 



'*» 



so that 



ty(g) 

dq 2 



i%v-M<ff)$'? 



g 2 i- 






l- 



4 g 2 



P + 



X 



2 \ 



'D' 



32 



er+A/ \g' 



■ «'(2Yp. 



32" \g : 



This gives on integration, since D/g 3 is independent of q, 



V(l-^M)dg» 



3g : 



*(1 



g 4 ) 



i 



_£ D 

4 g 2 



P + 



2 



9 



32\ cr + A, 






v \2 



ft 1 -*)®)'* 



• * 



• • 



• (121) 



The value of q 2 on the right of this equation is the root of equation (117), <j>(q) 
being given by equation (118) in which/ is in turn given by equation (119). Thus 
q 2 will be obtained by the elimination of <p (q) and f from equations (117) to (119), 
omitting terms in e 2 . From these three equations we obtain 

q 2 f = ~i> (q) = — [P-i/T>P + ^/ 2 D 2 P], 
so that / is a small quantity of the order of e, and, omitting terms in e 2 , 

<P (q) = «P 

which is independent of q. Equation (117) now gives, as the value of q\ 



g 2 = 2 






a 2 + \ 



+ eP. 



• • 



• * 



. (122) 



The coefficient of e on the right hand of equation (121) now reduces to 



L 2Lj 



X' 



(% ~1~A. _ 



1-2 



X 



2 \2 



a 2j r\ 



1 ■*-* 

*q 2 



X' 






\2 



64 



,2 



of + \J \q' 



9 6 






X 



2 \3" 



er + \/ J\g 



D 



,2 



P 



The value of P is given by equation (56), and that of D/q 2 by equation (120). 
From these we obtain 



^P = 2(- 2 - -±- )(l2Lf-W + 4mr + 4p), 
q 2 \a 2 a 2j r\J 



,q 2 



\2 



J> ___ V 

_J_ ■■ ■ ■ — ^y 



24L — - 



a 



2 



1 \ 2 7/I 

' +16Z i 



1 



2i B Z 
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whence we readily obtain, writing A for a 2 + X and so on, 



1 00 



AE,= 



J\-/ \ 



9 9 






wabc \ \ ( 2 ~ ) ( 2L ^7+ 22? J—- + 22p — 5 j 



XX 



B 2 C ; 



A 2 / 



H-g- 



l-IZ* 



9\ 9 






s <M 12L A 2+4 4 2+4 ' m $ + 4 ^ 



96 






61"( 2 — ) + X92 ( 2 



XX 



2\3 



242-^-- +162 A 



« 4 A : 



&VBC 



d\ 

A 



. (123) 



31. From § 8, the value of V,- (l) is given by 



/•OO 

V,(l) + eAV f (l) + e 2 aV,(l)=-,ra6c [/+^(l)] 

Jo 



dX 

4mmm 



• • « 



(124 



in which the value of $ (l) is now (cf. equations (12) and (47)) 



!- r / 3 D 4 P 2 



*(l) = e[P-i/DP +1 fc/ a D 8 P] 

-ie 2 [DP 2 -|/D 2 P 2 + r k/W- ,, - 15 

+ e 2 [Q-|/DQ+^; /2 D 2 Q-Wo4/ 3 D 3 Q + ...] 

+ terms in e 3 , &c 



* * 



. (125) 



Thus the value of SY { (l), being the coefficient of e 2 on the right hand of equation 
(124), is given by 

dX 



Jo 






Jo 



DP 2 -i/D 2 P 2 + t l/ 2 D 3 P : 



9 



isWW 



A 

i dX 
A 



32. Collecting and rearranging terms equation (115) now becomes 







[Q-4/DQ+ A/ 2 D 2 Q-^4/ 3 D 3 Q+ •••] f -^Qo 

= 0(y-2) [i(y-8)(z gV-[F + ( r -2) F'J (z£ 



+ i 
8 



,oo 



[DP 2 -i/D 2 P 2 + T i ¥ / 2 D 3 P 2 -Wx6/ 3 D 4 P 3 ] 





A 



AE ?: 
wahc 



+ Sn (x 2j rt/ 2 ) 



• * • * 



(126) 



This is an equation for the determination of Q, all the terms involving Q being on 
the left hand, and all the terms independent of Q on the right. The terms on the 
right are all of degrees 6, 4, and 2 in x, y, z, so that the value of Q will clearly 
consist of terms of degrees 6, 4, and 2 in x 9 y 9 z. Further, from the linearity of the 
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equation, it is clear that the value of Q will be the sum of a number of contributions 
arising from the separate terms on the right. 

Consider now what is the contribution from the first term on the right. Or, in 
other words, consider what would be the solution if the whole equation were 
reduced to 

'[Q-i/DQ+^/ 2 D s Q-^W 3 D 3 Q+...]^-eQo 

LI 

= i0(y-2)( y -3)(2^' 3 
The solution would be 

Q =-i( r -2)(y-3)(2^Y, (127) 



for, in accordance with the principles arrived at in § 14, this value makes the integral 
on the left vanish, and the equation is then satisfied. 

Similarly the contribution from the term in P" on the right is 

Q = ( r -2) 2 P // ^^), (128) 

for P /; is equal to ^ (2x 2 fa 2 ) 2 > and so the integral on the left again vanishes. 

The contribution from the term in P' in. equation (126) cannot be so simply 
evaluated. If, however, we are content to use Approximation B for P 7 , then P' 
becomes proportional to P /; , and an approximate solution is 

Qo = (y-2)F (2^. ........ (129) 

On referring back to § 16, it becomes clear that the accuracy of the approximate 
solution is of the same order as that which we previously called Approximation A. 

Again, if we use Approximation B for P , the contribution from the whole second 
line in equation (126) will be 

for the whole second line in question now represents merely the second order terms in 
the potential of a thick homeoidal shell, and so vanishes. 

We proceed next to the contribution from the term in AE* From equation (123) 
it appears that if we use Approximation B for the values of L, M, N, ..., and also use 
equations of the type of (69), the first line in AE t - ? as given by equation (123), 
becomes a function of (2x 2 /a 2 ), while the remaining lines vanish. 

It follows that an approximate solution is 

4o = ~^atee' 
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There remains the contribution from the term Sn(x 2 +y 2 ) on the right of equation 
(126). This contribution will consist of second degree terms in Q , which represent 
merely a step along the fundamental series of ellipsoids to allow for the altered 
rotation, and the requisite terms can be easily calculated. 

Collecting all the approximate partial solutions which have been obtained, we find 
the complete approximate solution 



Q = -(y-2) 



i (y-8) (Z £)' - [F + (y-2) F'J (2 £ 



H 7-^+ terms in Sn (130) 

Clearly this approximation is of a degree of accuracy comparable with that of our 
previous Approximation A for P , with which it may be compared. This former 
approximation can be put in the form 

P = |- ( y »2) (2 —}\ 2 + -i- + terms in An (131) 

33. A less good approximation, comparable with the former Approximation B, can 
be obtained by further simplifying equation (130) by the help of approximate 
equations such as (69). The value of Q simplified in this way is found to reduce to 
a function of (Zx 2 /a 2 ), so that the whole solution becomes ellipsoidal. 

34. The accurate solution for Q may be supposed to be 

Qo = Q'o+Q"o(y-2), 

where Q"o(y — 2) is the part of Q which is accurately given hy formulas (127) and 
(128); thus 

Q /, o = »(y-2)+«(2^' 

I have calculated the value of Q' accurately for one configuration only, namely, 
the ellipsoidal point of bifurcation. At this point a = b and x, y enter only through 

x 2 + y 2 . Let us write 

x 2 + y 2 = w 2 ; g 2 + £ 2 = m 2 . 

The general value of Q will be of the form 

QJA) a , OK5 4. o . ■ O Jl o a U a r a ZjS 9 9 . t A , ^U o ZjV <) 

a u a s c 4 a 4 c 8 c 12 a 8 a 4 c 4 c s a 4 c 4 

We may further put 

K = R' + K"(y-2), &c, 

then the values of B/ ; , &c. , are 

fi" = [i(y-2)+i]a 6 ,&o., 



OF ROTATING COMPRESSIBLE MASSES. 



187 



and the value of Q' is 



M 



T>/ OCJ/ 

it « . Okj 



a 






10° + 



a 8 c* 



^V+.., 



so that 



Q ; = RV + 3S'wT' + 3TW+UT J + ^ . (132) 



Let us now assume that 



,00 



(Q'-i/DQ'+ hP D 8 Q'-W W / 8 D 3 Q') 



,\ c?x 



A 



- h 1 (5w (i -90w i z 2 + 120w 2 z 2 -16z i )+h,{Hw i -24 : w 2 z 2 +8z i )+h R (w 2 -2z 2 ), . (-133) 

this being the most general form possible, since the integral in question is necessarily 
a zonal harmonic. In calculating the value of DQ', we may use the transformation 



r> x x 3 („ d U x a " 

a 2 A © 3rar \ Star/ c 2 C 3£ 2? 



and on equating coefficients in equation (133),. we obtain 



,oc t— 



5fh = 






A« 4 A 6 L 2 A 3 6R/ + c 2 C 6S / + 64 A* VA 2 



XS ,576R'+ — ^192S'+-£;72T' 



Uy C x\.vy 



4i^2 



c*U 



1 



ZdU4A \C£ xx 



fV-.2304R' + ° A ° 



aVA a C 



84S' + -J^U 288T' + ~ 720U' ) 



a¥AC 2 



c 6 C 



dX 



» * » 



(134) 



Olb<), — 



/k 



»oo 



r ' * / X 16r , + 4.4s'U4-/4r36R / +-46S / 



A 4 A \€(/ A 



c 2 C / 4A 4 VA 



c 3 C 



+ ^i-» i'-^;64r / + ?}' „■ 1 6s' + -^24t' 



6 4 A 2 VA 2 

1 / X 2 ^ 
32A 3 VA 2 

3 



576R'+ 



aVAC 

2X 2 



c a O 



X 5 



_96S' + -_^_72T / 

ffVAO e 4 C 



+ I X3 2304R' + 3X3 

+ 2304A 2 WA 8 + 



3X : 



X r 



4^2 A2p ^^ * ^2^4 



a c 



d\ 



384S' + -^^288T / +'~720UM ^ 



,00 



2u' 

A 2 4 A WA 



1 / X 8u'+-^4v') + -_L( -^16^+^48' 



c 2 G ^ / ' 4 A 2 Va 2 A 



./ 



32AWA ab4 aVAC H c 4 C a 



H —rrrl , , „ 576R ; H — o o a ; q96fe> H — r^r72l 

64A 2 VA 2 a 2 c 2 AC c*C 2 

X ^ r 2304R / +^SU384S / + 3X3 



"ZoU4A \Cfr A 



4.2 A 2 



a c 



a VAC 8 



288T / +-^r720U / 



c 



6 C 3 



* * # 



c?x 

. (136) 



188 



MR. J. H. JEANS ON THE CONFIGURATIONS 



35. The value of P is 



whence we find 



P = LTS* + 2mwV + N^ + 2pm 2 + 2rt; 



2 



DP 3 = EwH3£wt 2 + 3&TffYH©r + rw 4 + 2SwT + tf + 2ttw 2 + 2bf, 



in which 



3a = ^A : 


x 64L 2 + -^ 5 


< 8Lm 


<S — 


48Lm 


4(4m 2 + 2LN) 


© = 


J/(4m 2 + 2LN) 


40mN 


Wi = 


1 6mN 


56N 2 


x = 


144pL 


2 (8pm + 4rL) 


s = 


8 (8pra + 4rL) 


6(4jpN + 8rm) 


t = 


4 (4pN + 8rm) 


120Nr 


u = 


32p 2 


8pr 


b - 


16jpr 


24r 2 



We may now assume (cf. equation (133)) 



(DP-i/D 2 F+ T h/ 2 D 8 F- g A 6 / 3 D 4 F) 

Jo 



i\ c?X 



y x (5w 8 - 90wV+ 120w¥- 16z 6 ) +j' 3 (3w 4 -24wV + 82*) +/ 3 (w 2 -23 2 ), . (137) 



in which, on comparing coefficients, we obtain 

% 1 / X „**» . X „«\ . 1 / X' 



5? 



U-4-36E+-^6sU * 



A 6 8A 6 V« 2 A 

1 / X s 



c 2 C 



192A"\a 4 A 2 
3X 8 



9216A 8 \o 6 A 8 



576E + 



3X : 



VAC 



« 



192$+-^-„72C 



X ; 



^r2304E+ J 2 " -384g + ^ 8 " 4 " ^„ 288& + -^720 



arc 



a 



»c 4 AC ! 



c 



8Q8 



• ■ 



c 4 C 2 

dA 

A 

. (138) 



There are corresponding equations for 3j 2 and j s> but these are not written down 
as they can readily be obtained from equations (135) and (136). To change 3h 2 into 
3/ 2? change all accented coefficients into black letter type, and change X, X 2 , X 3 
wherever they occur explicitly, into ^X, |-X 2 , 3-X 3 respectively. The same procedure 
changes h% into j\. 

The values of these coefficients at the ellipsoidal point of bifurcation can be 
calculated from the material already provided in § 18 and the values of p and r 
given in formulae (113). Using the exact solution throughout, I find 

5j\ - 0*01931 (y-2) -0*01473, 

3/ a = 0*01080 (y-2) +0*00576, 

i 3 = 0*5715 (y-2)+0'4094, 

all terms in (y- 2) 2 vanishing in accordance with the result of § 14. 
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36. We turn now to the evaluation of AE^. The potential of the whole mass at 
an internal point, Y { , is given, as in equation (20), by 



Vi = /t>o(V t .(l)+eAy < (l)+ ...)-U-<r)(E i + eAE i + ...). 



We have also 



V 2 V,- -iTTp 



4tt/0 + 4 7T (/o — cr) ( 2— j +eP + e Q ~f ... 



whence, by comparison 

V»(AE,.) = -4,rP t 
We may assume 



-4tt 3 it 4- -r-y^ 4- - T r+--^w 4-3-2 
\a 8 a 4 c 4 c 8 a 4 c 4 



. (139) 



AE,- = - 4-7T ftw 6 4- 3k 2 w' k z 2 4 3k s wV 4- & 4 2 6 4- & & w 4 4- 2k 6 w 2 z 2 + k 7 z { + 2& 8 w 3 4- 2^ 



in which equation (139) necessitates the following relations between the coefficients 



r 

36^ + 6* 



2 ™ 8 

or 



< 



48L4-36L 



2m 






? 



12&. + 30& 



N 



*- c 8 



/" 



16& 5 + 4& 6 = ^?, 



^ 



8* 6 4-1 "2ife 7 



2r 



V 



iu/{Jg 4* AJq — U. 



By comparison with equation (123) 



ki = Ja6c 



LJ 



2L 



m 



A 5 



a 



2" X A 5 



2c 



L 
a 4 



m 



2 1 a 4 c + x 9 2 ( ^ 4 3 H A 5 4-32 "33 H A 4 C 4- 2 4 -j H A s c2 



a c 



N 

4 



c 



/6*k 



iabc 



T P T 1 ^* T 

2p J A 3 g 1-A 3 ~" 2" ~2 1 A 2 C 

a C 



2k 



8 — 4 



x a6c 



2Aj 



a' 



m 



L 



m 



N 



i~ Ia 3 + 1 ~3 Ia 2 c ~~* 6 4 ( 64 3 H A s 4-32 33 H A 2 C 4-243 H AC2 



a' 



2 2 
arc 



By direct computation from these values 1 find 



ik x 



0*019734 (y-2) -0'021633, 
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abc 


j \ji.vt o*± yy~ 


4y^ 5 ___ 
abc 


-0 - 007359, 


8k 8 _ 
abc 


-0'017022. 




2 c 
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37. We now return to equation (126). The solution has been assumed to be 

Qo = Q ; o+(y-2)Q" , 

in which (y — 2) Q /; is the contribution from the terms on the right 



0(y-2) 



X 



2\Z 



x d 



l( y „ 3 ) X ± _( 2 ) F / o s* 



a 



a 



It follows that Q' must satisfy the equation 



>oo 



f 1 

~~4 



/DQ'+^/ 8 D ! 'Q'-^ ¥ / s D 8 (y+ ...] ^-0Q', 



X' 



a 



.00 



(y-2)F 2± ) + i [DPM-/D 2 P 2 +.-.] 



dX 



^0 



s. 



irahc 



+ Sn(x 2 + y z ). 



Substituting the values which have been assumed for the various terms in this 
equation, it becomes 

M5w 6 -90wV+120wV-16^) + M 



e 






2„2 i t 4 2U 2 , ^V 2 



War +-2M — 7- tr + 
c 8 a 4 



z 



(y — 2) ( — + ~ ) ( — W L + 



a 2 c 2 / \<r 






+ i D'i (5w 6 --90wV + 120wV- 16s 6 ) + j 2 (3w 4 ~24wV + 8z 4 ) + j 3 (iv 2 -2z 2 ) 



4 



abc 



\k x tv Q + 3k 2 w 4 z 2 + 3& 8 wV + \z* + & 5 w 4 + 2k 6 w 2 z 2 + & 7 2 4 + 2& 8 w 2 + 2& 9 2 ~ 



_u 



<5n(aj 2 + 2/ 2 ). 



8 8 « 



(140) 



Equating coefficients of w 6 , wV, wV and £ 6 , we find 



5^-0 



/ -(y~2)L / a ! 



a 



12 






abc J 



9 • • 



. (141) 



90h!- e 



3S / -(y-2)(2wV + LVy 



a 8 c* 



X LiKv 



*jLj 1+ ±^, . . . (142) 



abc 



120 hy- 8 



'3T'-(y-2)(NV+2mV)" 



aV 



120^ + ^ 



• • » 



(143) 



16^-0 



V ~(y-2) Wc 2 )' 



a 



12 



¥-y 1 + i l t - • • • (144) 



abc 
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On multiplying equations (141) and (142) by 18, 1, and adding, we obtain, with 
the help of the relations of § 36, 



18 



B/~( y -2)Lyi ^ rsS^Xr^H^mV + LV)" 



a 



12 






aV 



2 L 



da 



,,8 



. . . (145) 



Equations (142) to (144), treated in a similar manner, give 

3S / -(y-2)(2mV + L / c 2 )n , „ r3T , -(v-2) (NV + 2mV) 



e*V 



+ 3 



4 8 

a <r 



2 m 



6a 4 c 



4^4 > 



(146) 



2 



'3T / -(y-2)(NV+2m / c 2 y 



aV 



+ 15 



'^-(y-^NV 



c 



12 



2N 



c 



8 ' 



. (147) 



These equations, taken with one other, suffice to determine B/, 8', T 7 , U 7 . For 
additional equation we shall use equation (141). 

The values of j\ and h x have already been given. For 5h x I find, by direct 
computation from formula (134), 

5L = 0'00152R / -0 o 01345S / + 0'03926T / --0-03813U / ? 



so that equation (141) reduces to 



-B/-(y-2)LV 



a 



12 



H-0*00152R / -0 , 013458 / + 0'03926T / -~0 , 03183U / 



0*02214 (y-2)-0-02347. 



• * 



Solving the system of four equations (145) to (148), I find 



E'-(y-2) LV _ A . 



a 



12 



0-04794-0'04309 (y-2). 



3S--( y -2)(2mV + NV) = . 1894 _ . 2573 ( 2)j 



a 8 c 4 



3r-( y -2)(NV+2mV) ^ . 4388 _ . 6707 ( 2 ) s 



a 4 c 8 



U'-(y 2)NV = '2605-0-5077( r -2) ; 



leading to the values 

E' 

3S' 
3T' 

U' 



0-4155 
0-1894 
0'0506 
0-00346 



-17799 (y-2), 
-1-4585 (y-2), 
-0-4124 (y-2), 
-0-0375 (y-2). 

A \J A 



A 



v 



(exact solution). 



J 



(148 



• » • * 



(149) 
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This may be compared with the approximation expressed by equation (130), which 
gives 



K' = (y-2)LV-4-&i a 



12 



ahc 



38' = (y-2) (2mV + L'c 2 ) - ~ U 2 



4 o7 . oV 




3T' 



U' = 



a&c 



(y-2)(NV+2m'c 2 )- 4-34 aV 



aoc 







12 



(y-2)NV--i-*A 

aoc y 



0-3981 -17691 (y-2), 
0-2013 -1-4659 (y-2), 
0-0452 -0-4090 (y-2), 
0-00355-0-0381 (y-2). 



>> 



s (Approxima- 
tion A) . (150) 



j 



The error of Approximation A is now something like 5 per cent, in the terms 
independent of (y — 2), and there is also an additional small error of less than 1 per 
cent, in the terms containing (y— 2). 

38. We pass now to the discussion of terms of degree 4. Equating coefficients of 
w* 9 w 2 z 2 , and £ 4 in equation (140) gives 



24h 9 -6 



3h 2 -0 


V -2 (y-2) pa 2 

« 8 


_ 3,y 2 + 4& 6 1 

8 ahc' 


~2s'-2 (y-2) (ra 2 +pc 2 Y 


_ 24/j ! 8& 6 
8 aftc ' 


a 4 c* 


8h 2 -6 


t'-2 (y-2) re 2 
L c 8 J 


_ sy 2 + 4* 7 

8 abc ' 



> 



* t 



(151) 



Following the method of the last section, the last two equations may be replaced 



by 



8 



r' -2 (y-2) pa* 2s'- 2 (y-2) (ra 2 +pc 2 )' 



a' 



+ 



4 4 

ere 



!£••■■ < 152 » 



'2s' -2 (y-2) (ra s +joc 8 )"l , „ ft'- (y-2) re 



4 4 



+ 3 



& 



2r 



(9c 



4* 



(153) 



The three equations (151) to (153) determine r', s 7 and t 7 . The values of 4k 5 /abc 
and of 3y 2 have already been given. For 3h 2 I obtain by direct computation from 
formula (135) 



OKI' 



3 A, = , 010972r / -0-064324s' + , 094925t 

+ 0-01031R / -0-007283S'-0-04422T' + 0-32343U', 

which, on inserting the exact solution (149) for R', S', T', U', becomes 



3h 2 = , 010972r / -0-064324s' + 0-094925t / + 0-00l79-0-00166(y-2). 
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The smallness of the last two terms of course measures the closeness of the 
approximation (130). 

Equation (151) now reduces to 



e 



"r' -2 (y-2) pa/ 



a s 



+ , 010972r / -0*064324s' + 0*094925t / 

= -0*00843 + 0*00310 (y-2). . . (154) 
On solving equations (152) to (154) directly, I find 

r> -2 (y-2) pa* = o- 2074-(r01055 (y-2), 
2 B'-2(y-2)(m'+,p ») = _ -07096 + 0-08440( y -2) ; 



aV 



t ' 2( y ^ 2)rc 2 = _0-i4008-0'02813(y-2) ! 

is 



whence the values of r', s', t' are 



? ^ ? 



r' = 0-08755 -010962 (y-2), 

s' = -0*01727 +0*04871 (y-2), }>(exact solution). . . . . (155) 

t' = -0*007862 + 0*02511 (y-2). , 

This may be compared with Approximation A, namely, 

which leads to the approximate values 

r' - 0*06590 -0"06509 (y-2)," 

s' = -0*00728 +0*02815 (y-2), ^(Approximation A), 
t' = -0*007929 + 0*02669 (y-2). 



The percentage error is large, although the absolute error in the coefficients is 
fairly small. This will be readily understood on noticing that under Approximation B, 
r', s' and t' would vanish altogether. 

39. Finally equating coefficients of w 2 and z* in equation (140) and making use of 
the relation k 9 = — 2k 8 , we obtain 

h-^4e = iJ a +^+Sn, (156) 

cr abc 

-2h s -^'-e=-h\-^ 057) 
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Multiplying by 2, 1 and adding, 



,/ ,„/• 



3n = -2d(^- + -A, . (158) 

an equation which may be used to replace (156). 
By direct computation from equation (136), I find 

h s = 0*19769u'-0*58215v' 

+ 0*06706r'-}-0*01164s / -0*61497t / 
+ 0-04270R / + 0'03778S / -0*07530T / -ri894U / 
- 0*19769u / -0*58215v / + 0*02348-0*11113 (y-2), 

so that equation (157) reduces to 

0*39538u' + 2*8142v / - 0*0212 + 0*0794 (y-2) (159) 

Equations (158) and (159) are adequate to give u' and v ; in terms of any assigned 
value of Sn. We require especially the solution at the actual point of bifurcation of 
the compressible series, and to obtain this the two equations must be combined 
with a third equation expressing the condition for a point of bifurcation. 

The third equation is 

S(a Q cJ AA )=0, 

or, transforming by the method of § 27, 

6u' T v ; T 2Ap 2 tt , 6pr yy 3r 2 rr 

a c a a c 2c 

Inserting numerical values, this becomes 

0'24919u' + 0'2Q878v / =: 0*00024. . ...... (160) 



From this and equation (159), I found by direct solution 

u' --0*00550-0*02651 (y-2) 
v' = 0*00778+0*03195 (y-2) 




161) 



and equation (158) now gives 

ck --0*01292-0*05495 (y-2) (162) 

The solution which has just been completed, combined with the two exact 
solutions (127) and (128) previously known, will give the exact solution for a 
compressible mass at the point of bifurcation at which the pseudo-spheroidal series 
gives up its stability to the pseudo-ellipsoidal series. 
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40. The rotation at this point of bifurcation is given by 



CO 



2tt/3 



n + e An + e* §n 

018712-0*04400 (^^)-[0'01292 + 0'05495 (y~2)] (^^Y. (163) 

\ Po / ~ V po 1 



The first point of interest about this equation is that the term in (p — a) is 
independent of y. This must necessarily be the case, for we have seen that for an 
incompressible mass, regarded as a special case of the compressible mass that has 
been under discussion, y — 2 is infinite while p —<r vanishes, and the product of the 
two remains finite. Thus for an incompressible mass (163) reduces to 

»72<7rp = 0'18712, 

which is the true value, but if there has been a term in (y—2) multiplying (p — or), 
equation (163) would have led to a wrong value. 
When y = 2, equation (163) reduces to 

co 2 /2ttp - 018712~-0'04400e--0-01292e 2 --.... 

Although we cannot be perfectly sure, the series is almost certainly convergent 
right up to the limiting case of <r = or e = 1. In this case, it appears to converge 
to a limit of about w 2 /27rp = 0*120, but unfortunately it is impossible to evaluate the 
limit for other values of y. 

A more important question is the relation between w 2 and p at the point of 
bifurcation. 

We have evaluated p as far as e 3 , but for comparison with the value of w given by 
equation (163), it will be enough to evaluate p as far as e 2 . Furthermore, to avoid 
a very complicated integration, we shall use Approximation B for the terms in e 2 , and 
so take 



(po — <r) 



Po — ypo — v 



x 
a 



, y , z , x " y J 

5 + - 6 T + ?+«»l^+y+-» 



^ 



C 



where y\ •==■ \{y— 2*75), this being the approximation given by equation (84). By 
a simple projective transformation, it follows that the mean density is the same as 
that in a sphere in which the law of density is 



■V 



A 



Po — (po~^)(— 2 + ej 7— ) • 



As far as terms in e, the radius of this sphere is given by r = (l—^erj)a, and the 
mass taken as far as e 2 is 



4 



7T 



'ro 



pr 2 dr 



■IT 



r< 



Po 



i \ i v 



Tr 



hw 



7a 



2 
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whence the mean density, as far as e 2 , is found to be 
From these equations we obtain 



, (164) 



O) 



2irp 



= 018712 + 0'06827e + (0'01602 + 0*07098(y-2))e 2 + (165) 



Ratio of Centrifugal Force to Gravity. 

41. [§41 revised August 31, 1918]. — -The ratio of centrifugal force to gravity is 
interesting only when it becomes equal to unity. When this occurs a stream of 
matter is ejected at the points at which centrifugal force equals gravity. Centrifugal 
force will first become equal to gravity at points which are furthest from the axis of 
rotation, and these will be points on the equator of the rotating mass. Considering 
a point on the axis of x> the condition for centrifugal force to be equal to gravity is 



or 



dYjdx + wX = 0, 
dQ/dx - 0. 



Since, by equation (30), = — ^~~ p y ~ l + a constant, this may equally well be 

1 



r 



expressed in the form 



dp/dx = 0. 



The general value of p is given by equation (45). On the ^-axis this reduces to 



1 — ev )j 



where 



F = 



/y» fc 

a" 



+ e 



L^ 4 2y>ff s 



<x 



a 



H-e 



2 



Ra? 6 . rx k . 2x\x' 



a 



12 



' + 



+ 



a c 



a 



~t~ • » • . 



(166) 



The intercepts on the axis of x are determined by the condition p = <r, and so are 
;iven by F = 1. The solution of this equation is found to be 



x_ 
a' 



= 1 



X , 2p 

p I — dL_ 



R _r 2u 

6 ' 4 ' 2 

a 3 a a 



,a 4 - a 3 / \ a 4 ' a 2 



4- 

I • • • » 



(167) 



The points on the #-axis at which cp/dx = are given by dF/dx = 0, so that it 
appears that dpjdx will just vanish at the extremities of the x-axis if d¥/dx vanishes 
for the value of x given in equation (167). 

The equation dF/dx = becomes, on division by 2x/a 2 



1-f 



V W7 a 2 



+ e 



2 



§R /V \ 2r /a; 

a 6 W j 



a 4 W/ 



2u 



+ ... = 0. 
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and, on inserting the value of x 2 fa 2 given by equation (167), this becomes 



1 + 



2L 2p 

a a . 



+ e 



3E , 2r , 2u 2L /L , 2p 

a 6 a 4 a 2 a 4 V a 2 



+ ... = 0. 



42. Let us examine in particular the special form assumed by this equation for the 
special configuration at which the pseudo-spheroidal form becomes unstable, giving 
place to a pseudo-ellipsoidal form. * Inserting the numerical values for L, p, R, &c, 
obtained in §§ 18-39, the equation becomes 

l+e[(y-2)-r0509] + e 2 [-Hy-2) 2 -0'4063 (y-2) -0'0510] + ... = 0. . (168) 

For a given value of y this equation determines the value of e for which centri- 
fugal force just outbalances gravity as the pseudo-spheroidal form gives place to the 
pseudo-ellipsoidal. 

For instance, for the value y = 2, the equation becomes 

l-r0509e-0'0510e 2 --.... 

The values of e obtained by using terms as far as e and e 2 respectively are 0'95 16 
and 0'9112. The true root is perhaps somewhere near e = T V Thus a mass of 
rotating matter obeying Laplace's law (y = 2) will throw off matter from its 
equator before reaching the ellipsoidal point of bifurcation if e > T 9 o, or if o- is less 
than x 1 ^/). 

The root e = T % agrees well with the corresponding quantity in the two-dimensional 
problem. For cylindrical masses obeying Laplace's law (y = 2) the problem can be 
solved exactly and the root is found to be e = 1, giving a- = 0. # 

Equation (168) may more usefully be regarded as giving y when the value of e is 
assigned. The only value of e which is of any astronomical interest is the value 
e = 1, the value for a mass whose density reduces to zero at the boundary. Putting 
e = 1, we obtain an equation giving a critical value of y. 

From terms as far as e only, the value of y is clearly enough 

y = 2'0509. 
On including terms as far as y 2 , this value becomes 

y = 21521. 

These values of y appear to be converging to a limit; we cannot state it with 
great accuracy, but we shall perhaps not be far wrong if we assume it to be y = 2*2. 

Here again we may compare the problem with the simpler two-dimensional one in 
which, as follows from what has already been said, the root of F = 0, when e = 1 , 
is y = 2 exactly. 

* 'Phil. Trans./ A, vol. 213, p. 471, 
VOL. OCXVIII. — A. 2 B 



198 MR. J. H. JEANS ON THE CONFIGURATIONS 

The Influence of other Physical Factors, 

43. The discussion from § 10 onwards has been limited to material in which the 
relation between pressure and density is of the form 

p = Kp y — cons., (169) 

an equation which implies that the mass is composed of matter of similar properties 
throughout. In an actual mass there must be at least some tendency for the 
substances of greater atomic or molecular weight to sink towards the centre, 
although it is far from clear to what extent this tendency will actually prevail 
under astronomical conditions. 

In equation (169), y is usually regarded as expressing the measure of a property 
of the material, being the ratio of the specific heats when the matter is gaseous. We 
have, however, only used this equation in its differential form 

dp = Kyp y ~~ l dp. . ......... (170) 

Now with any structure whatever of the mass of matter under consideration the 
surfaces of constant pressure and density must coincide. Our use of equation (170) 
has in no sense referred to a change of state of a single element of matter, it has 
referred to a passage from one surface of constant density to the next. Thus, in 
strictness, y refers to the arrangement of the matter and not to a physical property 
of the matter. In the case so far considered in which the matter has been supposed 
to have the same physical properties throughout, the two meanings become identical, 
but when the mass under discussion is an aggregate of imperfectly mixed types of 
matter, it becomes important to distinguish clearly between the two meanings. 

When p and p both vanish at the surface of the mass, equation (170) is equivalent to 

dlogp = -dlogp, (171) 

y 

where the differentials are supposed to refer to passage from one surface of constant 
density to the next. In passing from one surface to the next, p will increase partly 
on account of the increase of pressure, but partly also on account of the higher 
atomic weight of which the materials of the layer of higher density may be supposed 
composed. Thus we may suppose equation (l7l) replaced by 

dlogp = (_ + -!) d log _p, (172) 

where y M is the value of y defined by equation (169) as a property of the material, 
while y A is such that the additional term on the right of equation (172) represents 
the increase in d log p owing to increase of atomic weight of the materials as we pass 
to layers of greater density. 
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In any arrangement of matter whatever, y A will have some value for each layer of 
constant density, being virtually defined by equation (172). There is no reason to 
expect that y A will be constant, or even approximately constant throughout the mass. 
But in the present complex problem we must be content to discover tendencies rather 
than exact results, and so may think of y A as a constant. An increase of atomic 
weight on passing to layers of greater density will give a positive value to y A , while 
similarly a diminution of atomic weight would give a negative value of y A . 

The whole of the foregoing analysis will now apply to the present case if we put 

1 ~~* . . . • * . . » • . • l J. i ij f 



It has been seen that for values of y less than a critical value, which may be taken 
(although only as a rough approximation) to be 2*2, it will be impossible for the 
rotating mass to assume the ellipsoidal form, and so impossible for it to separate into 
two detached masses. The condition for assuming the ellipsoidal form now becomes 

7m Ta ^ ^ 

Thus, if y A is positive (atomic weight increasing towards centre), the critical value 
of y M will be greater than 2*2 ; with y A negative, the critical value of y M will be less 
than 2*2. It is possible for a heterogeneous but perfectly incompressible mass 
(y M = go) to fail to attain the ellipsoidal shape if y A is small enough — i.e. 9 if the 
layers increase sufficiently rapidly in density from increase of atomic weight alone. 

44. In the foregoing discussion we have considered only the effects of continuous 
changes of atomic weight. Entirely confirmatory results may be obtained from a 
consideration of excessively abrupt changes. 

Let us consider only the simplest case in which the mass is formed of completely 
separated layers of two substances only, one very light and the other very heavy — 
to make the picture definite, let us think of an inner mass of iron vapour surrounded 
by an atmosphere of hydrogen, and let us suppose the density of the hydrogen may 
be neglected in comparison with that of the iron, so that the gravitational field will 
be the same as if the iron alone were present. 

For any given value of oo we can draw the equipotentials Q = cons., and the 
boundary of the mass of iron must be one of these, say = 0^ Proceeding outwards 
and drawing the exterior equipotentials we must in time come to one having a double 
point, say Q = C x , and at the double point on this equipotential centrifugal force will 
be exactly equal to gravity (dQ/dn — 0). The space between the equipotentials Q = G x 
and Q = C x can be filled with hydrogen without matter being thrown off by rotation, 
but clearly no greater volume of hydrogen than this can be retained. 

As the rotation go increases from zero upwards, the equipotential Q = C z starts 
from infinity and moves inwards, so that the maximum volume of hydrogen which 

2 D 2 
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can be retained continually decreases. With an allotted amount of hydrogen, the 
ellipsoidal point of bifurcation may or may not be reached before matter has begun 
to be ejected from the equator of the figure. 

The simplest case for accurate discussion occurs when the core is treated as 
incompressible. In fig. 3 the thick curve represents the cross-section through the 
axis of rotation of an incompressible mass at the Maclaurin-Jacobian point of 
bifurcation,* while the thin curve represents the equipotential Q = G z9 which has a 
series of double points round its equator. If the volume of light material 
surrounding the core is just equal to the volume between these surfaces, matter will 
begin to be thrown off at the equator precisely at the moment at which the 
ellipsoidal point of bifurcation is reached. If the volume of light matter is at first 




Fig, 3. 

greater than this critical volume, matter will be thrown off equatorially before the 
ellipsoidal point of bifurcation is reached, the amount ejected being such that the 
volume is just reduced to the critical volume when this point is reached. Conversely, 
of course, if the volume of light matter is initially less than this, the point of 
bifurcation will be reached before any matter is thrown off equatorially. 

A rough estimate shows that the critical volume, when the matter is incompressible, 
is about one-third of the volume of the core. For a compressible core for which 
y > 2*2, it is of course less, becoming equal to zero when y = 2*2 (approximately). 

This completes our collection of theoretical results. They may now be recapitulated 
and discussed with reference to actual astronomical conditions. 

Summary and Discussion of Results. 

45. Our discussion began with a general survey of the types of configurations 
which can be assumed by compressible astronomical masses in rotation. The result 
announced in a previous paper # that the incompressible mass provides a good model 
from which to study the behaviour of compressible masses has on the whole been fully 
confirmed. A mass of incompressible matter, shrinking while rotating, will assume 
first the shape of a spheroid, then that of an ellipsoid, afterwards becoming unstable 

* 'Phil. Trans.,' A, vol. 213, p. 457. 
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and probably dividing into two detached masses after passing through a series of 
pear-shaped configurations. In the present investigation it has been found that, 
except for an alternative possibility to be discussed later, a shrinking compressible 
mass will experience a very similar sequence of changes. It will pass through 
a sequence of figures which from their similarity to the spheroids and ellipsoids of 
incompressible masses, may be described as pseudo-spheroids and pseudo-ellipsoids. 
The pseudo-spheroids become unstable when a certain degree of flatness is reached 
and give place to a series of pseudo-ellipsoids, these in turn become unstable when a 
certain degree of elongation is reached and give place to a series of pear-shaped 
figures which probably end by fission into detached masses. 

46. It may simplify the presentation of the detailed results obtained, if we confine 
our attention at first to the innermost strata, this being supposed to mean strata so 
near the centre (which is taken as origin) that terms in x*, x 2 y 2 , &c, may be neglected 
in comparison with terms in x 2 , y 2 , z 2 . It has been found that these strata will be 
first truly spheroidal, then truly ellipsoidal, then pear-shaped, each change of figure 
corresponding, as in the incompressible mass, to a passage through a point of 
bifurcation of the system as a whole. 

The problem has been studied in detail for a pressure-density law of the form 

p = K py~ cons., (175) 

and special attention has been paid to the exact position of the configuration at which 
pseudo- spheroidal configurations give place to pseudo- ellipsoidal ones. This configu- 
ration we may refer to as the ellipsoidal point of bifurcation ; its special importance 
will appear later. The corresponding configuration for an incompressible mass is the 
Maclaurin-Jacobian point of bifurcation ; at this point the semi-axes a , c and the 
angular velocity w are given by 

a = 11972, c = 0*69766, co 2 /2tt = 0'18712/>, 
so that the equation of the boundary is 

x 2 +y 2 z 2 _ - 

2 ' 9 • 

For a compressible mass, the equations of the innermost strata of constant density 
at the ellipsoidal point of bifurcation have been found to be 

^^-[l~0-016307€-(0-00384 + 0-01850(y-2))€ 2 ~...] 

t 

+ ^ i [l + 0'056346 + (0-01598 + 0-06565(y~2))6 2 +...] = e ^ = ^, . (176) 

C l p — or 

where p is the density at the centre, a- that at the surface, and e stands for (p —cr)/pw 
For the important case of a gaseous mass, e is of course equal to unity. 
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Thus at the ellipsoidal point of bifurcation, the innermost strata of equal density 
are of flatter shape than those for an incompressible mass, showing that compressibility 
tends to postpone instability for the spheroidal mass. 

The value of w, the angular velocity at this point of bifurcation has been found to 
be given by 

c t ) 2 /27r / o = , 18712-0 , 04400e-(0-01292-h0 , 05495(y-2))€ 2 -... ? . . (177) 
or, if we evaluate co 2 in terms of p, the mean density of the whole mass, 



co 2 



/2tt/> = , 187l2 + , 068276 + (0-01602 + 0-07098(y-2))e 2 -..., . . (178) 



the coefficient of e 2 now being only approximate. We notice that afjiirp is greater 
for a compressible mass than for an incompressible mass, so that again compressibility 
may be said to postpone the instability of the spheroidal form. # 

Equation (176) applies only to the innermost strata for which a? 4 , &c, may be 
neglected. The equation of the outer strata, both at the point of bifurcation and 
elsewhere, are found to contain terms of degrees four and higher, there being terms of 
degrees four and two multiplied by e, terms of degrees six, four and two multiplied 
by e 2 , and so on. The terms of degrees six and four have been calculated for the 
ellipsoidal point of bifurcation and the terms of degree four for the pear-shaped 
point of bifurcation. 

The presence of these terms destroys the spheroidal or ellipsoidal shape of the 
outer strata. In general it is found that the outer strata are more lens-shaped than 
the inner strata when these latter are spheroidal, and more spindle-shaped than the 
inner strata when these latter are ellipsoidal. The lens- shaped form of the outer 
strata may go so far that the outer boundary develops a sharp edge. When this 
occurs, centrifugal force is exactly equal to gravity at points on the periphery of the 
lens, and any further increase in the rotation of the mass results in matter being 
thrown off from round this periphery. Similarly the spindle-shaped figure may 
develop sharp ends, in which case matter will be thrown off here also. 

47. Consider now a gradually shrinking mass of gas or other compressible matter, 
the rotation increasing as the shrinkage proceeds. For a very slow rotation the 
strata and the boundary will all be spheroidal. As the rotation increases, the 
boundary departs more and more from the spheroidal form, taking a series of forms 

* All terms in a> 2 /27rp are positive because y- 2 is necessarily positive; for, as we shall see, if y<2, the 
compressibility so far postpones the occurrence of the ellipsoidal point of bifurcation that it does not occur 
at all. Incidentally equation (178) has an important bearing on the origin of the solar system. It shows 
that for every mass which has broken up by fission (o 3 /27rp must at some time have exceeded the value 
0*18712. This provides the keystone, which has so far been wanting, in an* argument I have given 
elsewhere ( l M. N. Royal Ast. Soc./ 77, p 191) to show that the solar system is very unlikely to have 
broken up by rotation alone. 



OF ROTATING COMPRESSIBLE MASSES. 203 

which we have called pseudo-spheroidal, these being more lens-shaped than true 
spheroids. There are two alternatives. It may be that the point of bifurcation will 
be reached before a sharp edge forms at the equator of the pseudo-spheroids, and 
if this happens the innermost strata will become ellipsoidal, and the outer strata 
and the boundary will be pseudo-ellipsoidal ; the mass will proceed towards the pear- 
shaped form and ultimate fission. But, instead of this happening, it may be that a 
sharp edge will be formed before the point of bifurcation is reached, and the mass 
will disintegrate through equatorial loss of matter. It is of the utmost importance 
to determine which of these events will happen first for a particular mass. 

48. Some information may be obtained from a general survey of the problem. 
The spherical solutions for a mass of matter at rest obeying the law p = Kp y have 
been investigated by Hitter, Darwin, Emden # and others. Excepting the special 
case of y = 2, in which the equation reduces to a linear equation, the solution can be 
expressed in finite terms in one case only, namely y = 1^. In this case the solution, 
first given by Schuster, t is 

P = A(l+r 2 /a 2 )~ l (179) 

This happens also to give the lowest value of y for which the mass is finite. For 
values of y less than 1-J-, the matter extends to infinity and the total mass is infinite ; 
when y = 1 ■£-, the matter extends to infinity but the total mass is finite ; when 
y > 1-J-, the matter is of finite extent and of finite total mass. As y increases the 
variations in density becomes less rapid and finally the value y = oo corresponds to 
an incompressible mass in which the density is uniform throughout. 

Now clearly a mass for which y = 1-J- will lose matter equatorially with even the 
slightest amount of rotation. For all except an infinitesimal fraction of the whole 
mass is concentrated in regions near the centre, and the potential near the edge may 
accordingly be taken to be M/r. The value of Q is therefore 

r 

The problem is now seen to be identical with one which has been studied by 
Roche. J There is a critical equipQtential on which double points occur all round the 
equator, and this is the equipotential 

Q = f (Mco) 1 . 

* E. Emden, 'Gaskugeln' (Leipzig, 1907), where references to the work of previous investigators 
will be found. 

t ' British Assoc. Beport,' 1883, p. 428. 

J " Essai sur la constitution et Porigine du Systeme Solaire," ' Acad, de Montpellier. Section des 
Sciences,' VIII., p. 235. A more accessible account is to be found in PoiNCARE's 'Lec/ms sur les 
Hypotheses Cosmogoniques,' 2nd edition, p. 15. 
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The equatorial radius is M^ -1 . By a simple integration the volume of this critical 
equipotential is found to be 

« ft 17 A pnq 2 ft — 3 

32ttMco- 2 , a * 7 cos 2 sin c?0 = 32*-M*>~ 2 x 0*0225466. 
Jo 4 cos 2 0-1 

This volume is of course equal to M/jo where ~p is the mean density, from which 

we find 

w 2 /2irp = 0*360746 (180) 

Since p = 0, the critical value of a> is also zero, and the innermost strata differ 
only imperceptibly from spheres. Thus, when y = 1 J, equatorial break-up occurs as 
soon as the mass is set into rotation at all, and therefore long before there can be 
any question of the pseudo- ellipsoidal form being attained. 

At the other end of the scale (y = oo) comes the incompressible mass, for which 
the ellipsoidal form is attained long before there is any question of equatorial 
breaking up. 

Thus for some value of y, intermediate between y = 1 § and y = oo 9 there must be 
a crossing over from equatorial break-up to fission through pseudo-ellipsoidal and 
pear-shaped figures. For a mass for which y has this critical value, the point of 
bifurcation is reached and the pseudo-ellipsoidal form assumed at the very instant at 
which equatorial break-up is about to begin. 

49. A comparison with the corresponding two-dimensional problem is of interest 

at this stage. Here again a solution in finite terms is only possible for one value of 

y other than y = 2, and here again this value happens to be that one for which the 

total mass is first finite, while the matter extends to infinity. The value in question 

is y = 1 and the solution is 

P = A(l+r 2 /a 2 )-\ (181) 

which may be compared with the three-dimensional solution (179). Again it is clear 
that, somewhere between y = 1 and y = oo ? there must be a critical value of y at 
which a transition occurs from equatorial break-up to fission into two detached 
masses. The two-dimensional problem can, however, be fully solved, and the 
critical value is found to be y = 2 exactly. # 

Seeing that the two problems run fairly close together in all their essential 
features, we might suspect that the critical value in the three-dimensional problem 
would not be very far from 2. An alternative guess might be y = 2|, this corresponding 
more closely to the two-dimensional value y = 2, since y = 1 J in the three-dimensional 
problem has been found to correspond to y = 1 in the two-dimensional problem. 

50. In the present paper the critical value of y has been shown to be determined 
in the three-dimensional problem by the equation 

l + e[(y--2)-r0509]+6 2 [|-(y-2) 2 -0-4063(y-2)-0'0510]-i-... =0, . (182) 

* 'Phil. Trans./ A, vol. 213,. p. 471. 
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in which e must of course be put equal to unity for a gas in which <r = at the 
boundary. The series on the right is probably rapidly convergent, but sufficient 
terms have not been calculated for the value of y to be determined with accuracy. 
Neglecting all terms beyond those written down, the root of equation (182), when 
e = 1, is found to be y = 21521, but the remaining terms appear likely to increase 
this value somewhat, and we may perhaps take y = 2*2 as an approximate value. 
This, we may notice, is just half-way between the two guess-values considered in § 50. 

When the density at the boundary is not zero, the critical value of y is less than 
this ; for instance, a critical value y = 1 corresponds to a value of e equal to about 
three-quarters. 

An approximately accurate drawing of the critical figure, when e = 1, is shown in 
fig. 4, the inner curves being strata of constant density. 

51. We have considered the effect of heterogeneity in the structure of the matter, 
and have found that a sinking of the heavier elements to the centre of the mass will 
result in an increase in the critical value of y. There is no limit to the amount of 




Fig. 4. 

increase that can be produced in this way, although naturally the amount of increase 
depends on the extent to which the Hght and heavy elements are separated and on 
the ratio of their amounts. As an illustration, we considered the extreme case of a 
core of heavy incompressible material which we called iron, surrounded by an 
atmosphere of much lighter material which we called hydrogen, the two elements 
being supposed to be completely separated. If the volume of hydrogen was initially 
greater than about one-third of that of iron, the composite mass set into rotation will 
first disintegrate through equatorial loss of matter. 

A drawing of the critical figure has already been shown on fig. 3 (p. 200) ; this 
may be compared with fig. 4. 

52. There is, however, a very essential difference between the failure of a uniform 
mass to attain the ellipsoidal form and the corresponding failure of a heterogeneous 
mass. Briefly speaking the former is permanent, while the latter is transitory. 

Consider for simplicity a mass of perfectly uniform gas — say helium for which 
-set into rotation and continually shrinking by cooling. The value of y is 



= 1 



2 
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certainly below the critical value (about 2*2), so that as the rotation increases matter 
will be thrown off from the equator before the ellipsoidal shape can be reached. No 
matter how much material is ejected in this way, the central mass remains a mass of 
helium with y = if, and so can never attain the ellipsoidal form. The mass will 
disintegrate completely by loss of matter from its equator, and there can be no fission 
into separate masses — unless, of course, the helium ultimately so changes its character 
that the average value of y becomes greater than about 2*2. 

Contrast this with the behaviour of the mass of iron and hydrogen already 
considered. The mass will lose matter by equatorial ejection, but the matter lost 
will consist entirely of hydrogen, and so this process will continually diminish the 
ratio of hydrogen to iron, and as the process continues the mass will continually 
approximate to a mass of incompressible iron. As soon as the ratio of hydrogen to 
iron is reduced to about one- third, the pseudo- spheroidal form becomes unstable, and 
the mass will assume a pseudo-ellipsoidal form. 

53. We are now in a position to follow the changes in a mass of gas whose rotation 
continually increases through shrinkage. At first we may assume the gas to obey 
the ideal laws, so that y will be less than if, and, as the rotation increases, a stage 
will be reached at which matter is ejected from the equator. Some of this matter 
will perhaps fall back on to the rotating mass, but some must also pass to infinity, 
this latter representing a real loss of mass and of angular momentum to the rotating 
body. The loss must be at such a rate that the figure of the rotating body always 
remains a pseudo-spheroid with a sharp edge, and the velocity of rotation remains 
always exactly equal to the critical velocity corresponding to this critical figure. 

When y is equal to 1-J, the critical angular velocity is given by equation (180), so 
that w 2 = 0*36x27r£. For the greatest value of y for which the pseudo-spheroidal 
figure is possible (about y = 2*2), the value of the critical angular velocity is given 
by equation (178), and this seems to be converging to a value not far from 
co 2 = 0'36 x 27r/o. We may perhaps conjecture that for all masses of gas which are 
throwing off matter equatorially the value of w is nearly equal to 0'36 x 2xp. 

As the mass shrinks p becomes greater, so that rotation becomes more and more 
rapid. A stage is reached in time in which the ideal gas laws no longer hold, owing 
to the distance apart of the molecules having become comparable with their diameters. 
The value of y now increases beyond its value for an ideal gas, and values of y 
greater than if become possible. When a value of y is reached which is about 
equal to 2'2 if the mass of gas is perfectly mixed, but may be greater if the mixing is 
imperfect, the pseudo-spheroidal form becomes unstable, the mass assumes the pseudo- 
ellipsoidal form and the process probably ends by fission into two detached masses. 

. 54. The value y = 2*21 has been found by Koch for air at a pressure of 100 atmo- 
spheres and a temperature of — 79°'3 C, the corresponding density being 0*23. 
Partly from this observational material, and partly from general theoretical 
principles, we may anticipate that the value y = 2'2 will be attained at a density 
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of about one-quarter of that of the substance in its solid state, and for the permanent 
gases this may be taken to be about a quarter of the density of water. 

Thus a mass of gas will lose matter equatorially until it has shrunk to a density 
of about a quarter of that of water, after which it will elongate and divide up. 
Assuming the relation w = 0*36 x 2irp, the mean density p and the period of rotation 
in days (P), will be connected, so long as the mass is losing matter equatorially, by 

the relation 

p = 0'035 -~P 2 , 

the density p = i corresponding to a period of about 9 hours and p = |- to a period 
of 6-J- hours. When this stage is reached the process of elongation followed by 
fission begins. Assuming that when fission is complete we have two stars of 
approximately equal mass and mean densities p = J revolving round one another 
almost in contact, the period of this system would be about a day. 

55. The critical density which we have conjectured to be about one quarter is 
perhaps not far from that of the average B-type star. # Thus, subject to the 
assumptions on which we have been working, fission ought to begin at about B-type. 
This is in very close agreement with the results obtained by Campbell in his 
" Second Catalogue of Spectroscopic Binary Stars." t 

It is not, however, in agreement with the results obtained by ShapleyJ in his 
" Study of the Orbits of Eclipsing Binaries." Adopting Kussell's view of stellar 
evolution, our result would show that giant stars (except, possibly, of A- type) 
should be pseudo-spheroids; only B and dwarf stars could form binaries. Shapley 
discusses 93 systems; 88 are of B or dwarf type, but only 21 have densities greater 
than 0'316, and only 57 have densities greater than O'l. For one star, W Crucis, 
of which the orbit has been determined by Russell, § the density of the brighter 
component appears to be of the order of 0*000002. We are led to inquire under 
what physical conditions, different from those we have assumed, it can be possible 
for fission to occur while the density is still far below the value to which our 

analysis has led. 

56. Consider the simplest problem of a spherical mass at rest, the equations of 

equilibrium being 



dp _ 3V _ 47rGp 
dr dr r 2 



j pr 2 dr, 

Jo 



where G is the gravitation constant. If p is the pressure at the centre and R the 
radius, we obtain, on integrating from the boundary to the centre, 

= 4ttG ■£ pr 2 dr dr. 



Jo r* 



o 



* Eussell estimates tne density of a giant A-type star at 0*1 (' Nature/ 113, p. 282). 

t ' Lick Observatory Bulletin/ 1910. 

% * Princeton Observatory Contributions, No. 3 ' (1915). 

§ * Astrophysical Journal/ 36, p. 146. 

2 E 2 
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If for any reason the ratio of increase of pressure to density is about that 
corresponding to the critical value y = 2*2, the variation of density, except near the 
surface, will not be very great. An approximation which will be accurate as regards 
order of magnitude at least will be 

p = 47rG^ 2 ~ r 2 drdr = -^ Gp 2 K 2 , 

Jo T Jo 3 

or, since f wpR 8 = M (the mass of the body), 

Po = iG{^Jmp^5xlO~ 8 Wp.. ...... (183) 

Assuming the gas law 

p = 5_T P B, r (184) 

m 

where R is the gas constant, and B a multiplying factor introduced by deviations 
from Boyle's law, we find for the temperature at the centre 

To = H & = 5 x 10- 8 ^ M^ (pM, 

Xt-D p Q XtJJ 

in which p/p may, in the case of equilibrium with y equal to about 2 '2, be put equal 
to about 0*6, giving 

Xq ~ o X 1U Tyfy P 3 * ' ' ' * (185/ 

The energy of radiation at the centre is o-T 4 per unit volume, where o- is Stefan's 
constant of which Kurlbaum's value is 7*06 xlO" 15 , and the pressure of radiation, 
being one-third of this, is -g-a-T 4 in all directions. 

Denote the pressure of radiation by p R and the gas pressure by p Q . If j9 , the pressure 
at the centre may be regarded as arising mainly from gas pressure, we have 

& = fe^ 4 = 4xlO- 88 (^)M a > (186) 

Pq Pq \±t±>/ 

which is independent of jo. 

For a gas of molecular or atomic weight 32 the value of w/R is about 4x 10~ 7 , 
while the value of B will not differ greatly from unity until very high densities are 
reached. Thus equation (186) becomes approximately 

^ = 10-M* (187) 

Pq 

It will only be when this ratio is small that the radiation pressure will be negligible 
in comparison with gas pressure at the centre of the star, and this fixes a limit to M 
of the order of 10 31 , independently of the density of the star. 
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Most stars whose mass is known* have a mass comparable with that of our sun 
(M = 2x 10 33 ), and for these the ratio (187) is far from negligible. Thus it appears 
that for stars of mass comparable with our sun, and of molecular or atomic weight 
about 32, the pressure of radiation will not be negligible in comparison with gas 
pressure ;t the equations of equilibrium must be replaced by 

and our calculation fails from equation (184) onwards. 

57. If we had taken a molecular weight 2, instead of 32, the value of p R /p a would 
have been reduced by a factor (16)~ 4 , and we should have had approximately 

£s = 10- 68 M 2 , 
Pa 

a ratio which may be considered small when M is of the order of 10 33 . Now whether we 
consider that radiation pressure is comparable with gas pressure or not, the tempera- 
ture at the centre of stars such as we have considered is of the order of 10 7 degrees 
Centigrade, and at such temperatures it seems probable that matter would to a large 
extent be broken up into its constituent electrons and nuclei. For purposes of 
calculation of gas pressure each electron behaves like the molecule of a gas, and the 
effect of electronic disintegration is to reduce the effective molecular weight. It is 
readily seen that when electronic disintegration is complete, a limiting effective 
molecular weight of 2 is reached for all substances except hydrogen.^ 

Further, radiation pressure when it is appreciable may be treated as arising from 
molecules of molecular weight zero, and so the effective molecular weight may be 
still further decreased. 

* The stars whose masses are known are bright binaries — binaries because there is no means of determining 
mass except by the mutual action of two bodies on one another, and bright because the fainter binaries 
escape observation. The conclusion of this paper is that bright binaries are binaries of large mass, the mass 
being of order of magnitude greater than 10 31 gm. This is borne out by the observed masses of those binaries 
which are bright enough to have attracted attention, but there is nothing to show that there are not a 
great number of less bright binaries of smaller mass. It rather appears as if the few well-determined 
masses of binaries are not likely to give a good sample of the masses of all stars. It is perhaps 
significant that Russell's well-known diagram of absolute magnitudes (' Nature/ vol. 113, p. 252) shows 
a range of something like five magnitudes (ratio 100 to 1) for dwarf stars of similar spectral type, 
suggesting a range of masses enormously greater than that calculated from observations on binary stars. 

t This agrees with the result stated by Eddington (' M.N., Eoyal Ast. Soc./ vol. 77, p. 16), but as the 
question is of some importance for our present investigation, I have thought it worth giving a separate 
discussion freed from the special assumptions of Eddington's paper. Eddington's statement that 
radiation pressure is practically negligible for dwarf stars does not appear to be altogether confirmed by 
my equation (187). 

J Cf. Eddington, ' M.N., Eoyal Ast. Soc.,' vol. 77, p. 596. 
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Thus the effective molecular weight, regarded as being determined by equation 
(184), may faU to very small values as we pass to the interior of a star. In § 43 we 
considered the tendency for greater molecular weights to occur in the interior of the 
star, and found that it would delay the stage at which the pseudo-spheroidal form 
would become unstable. It now appears that, in order to represent actual conditions, 
we should have examined the reverse tendency. Effectively lighter molecular 
weights may be expected to occur in the interior of a star, and this will cause the 
process of fission to begin at lower densities. In a general way we may expect that 
stars of greatest mass will begin the processes of fission in the earliest stages of their 
careers. If so, such a star as W Crucis, assuming that its two components have 
been formed out of the fission of a single body, must be a star of very great mass 
indeed. 

58. To sum up, we have found that a star of small mass (say xVth that of our sun), 
in which there is no great amount of atomic disintegration, and in which pressure of 
radiation does not play a prominent part dynamically, will not begin to break up into 
a binary star until it reaches a density of from one-quarter to one-half that of water. 
In more massive stars there will be considerable atomic disintegration, and pressure 
of radiation will be dynamically important. Such stars will break up at lower 
densities than smaller stars. 



